A method of generating a four dimensional seismic signal based on multiple sets of seismic data representing a subterranean formation. The method can include generating a tomographic velocity model based on a first set of raw seismic data and determining at least one Green's function based on the tomographic velocity model. The method can include generating a first image of a target region based on the first set of raw seismic data and the at least one Green's function. The method can include generating a second image of the target region based on a second set of raw seismic data and the at least one Green's function. The first and the second images can be compared, and a four-dimensional seismic signal can be determined based on the comparison.

A method of generating a four dimensional seismic signal based on multiple sets of seismic data representing a subterranean formation. The method can include generating a tomographic velocity model based on a first set of raw seismic data and determining at least one Green's function based on the tomographic velocity model. The method can include generating a first image of a target region based on the first set of raw seismic data and the at least one Green's function. The method can include generating a second image of the target region based on a second set of raw seismic data and the at least one Green's function. The first and the second images can be compared, and a four-dimensional seismic signal can be determined based on the comparison.

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.

An actuation location for actuation of a first source coupled to a first marine survey vessel relative to a position of a second marine survey vessel towing a receiver to enhance illumination of a subsurface location can be determined based on a survey route of the second marine survey vessel and a priori data of the subsurface location. The first marine survey vessel can be navigated along a survey route of the first marine survey vessel to the actuation location during a marine survey by changing at least a cross-line position or an in-line position of the first marine survey vessel relative to the survey route of the second marine survey vessel.

The present invention is related to a specific de-migration method based on ray tracing algorithms characterized in that the interpolation procedure involved in the computation of the travel time required by the de-migration is being modified. The interpolation according to the invention obtains an accurate travel time for those rays departing from sources being interpolated.

Systems and methods for estimating orthorhombic anisotropy parameters of subsurface rock layers are provided. An initial three-dimensional (3D) model of layers in a subsurface formation is generated. Particular combinations of seismic source and receiver locations associated with a vertical seismic profile (VSP) survey of the formation are selected based on a simulation of seismic wave propagation through each layer of the initial 3D model. A global inversion is performed using data points selected from travel time data associated with the VSP survey in order to estimate different sets of anisotropy parameters for the layers of the initial 3D model, where the selected data points correspond to the selected source-receiver combinations. The initial 3D model is refined with an optimal set of anisotropy parameters selected from the estimated parameter sets. The refined 3D model is provided for seismic analysis and well planning operations to be performed for the subsurface formation.

Systems, computer readable, and methods concern receiving seismic data representing a subsurface volume. The method also includes determining, for the seismic data, analysis coordinates as a function of time. One or more of the analysis coordinates may vary in position over time. The method includes performing at least one of an interpolation or regularization process on the seismic data based at least partially on the analysis coordinates. The method also includes outputting a result of the at least one of the interpolation or regularization process.

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.

Systems and methods for generating a three-dimensional image of a proppant-filled hydraulically-induced fracture in a geologic formation are provided. The image may be generated by capturing electromagnetic fields generated or scattered by the proppant-filled fracture, removing dispersion and/or an attenuation effects from the captured electromagnetic fields, and generating the image based on the dispersion and/or attenuation corrected fields. Removing the dispersion and/or attenuation effects may include back propagating the captured electromagnetic fields in the time domain to a source location. The image may be generated based on locations at which the back propagated fields constructively interfere or may be generated based on a model of the fracture defined using the back propagated fields.

Aspects provide for methods that successfully evaluates multiple compressional and shear arrival events received by a sonic logging tool to evaluate the presence of structures, such as shoulder beds, in downhole environments. In particular, the methods described herein enable automated determination of properties of laminated reservoir formations by, for example, enabling the automated determination of arrival times and slownesses of multiple compressional and shear arrival events received by a sonic logging tool.