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.

A method is provided for identifying and characterizing structures of interest in a formation traversed by a wellbore, which involves obtaining waveform data associated with received acoustic signals as a function of measured depth in the wellbore. A set of arrival events and corresponding time picks is identified by automatic and/or manual methods that analyze the waveform data. A ray tracing inversion is carried out for each arrival event (and corresponding time pick) over a number of possible raypath types to determine i) two-dimensional reflector positions corresponding to the arrival event for the number of possible raypath types and ii) predicted inclination angles of the reflected wavefield for the number of possible raypath types. The waveform data associated with each time pick (and corresponding arrival event) is processed to determine a three-dimensional slowness-time coherence representations of the waveform data for the number of possible raypath types, which is evaluated to determine azimuth position and orientation of a corresponding reflector, and determine the ray path type of the reflected wavefield. The method outputs a three-dimensional position and/or orientation for at least one reflector, wherein the three-dimensional position of the reflector is based on the two-dimensional position of the reflector determined from the ray tracing inversion and the azimuth position of the reflector determined from the three-dimensional slowness-time coherence representation. 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.

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 and a system for implementing the method are disclosed wherein the seismic input data and land acquisition input data may be obtained from a non-flat surface, sometimes mild or foothill topography as well as the shot and receiver lines might not necessarily be straight, and often curve to avoid obstacles on the land surface. In particular, the method and system disclosed, decomposes the cross-spread data into sparse common spread beams, then maps those sparse beams into common-spread depth domain, in order to finally stack them to construct the subsurface depth images. The common spread beam migration and processing have higher signal to noise ratio, as well as faster turn-around processing time, for the cross-spread land acquisition over the common-shot or common offset beam migration/processing. The common spread beam migration method and system disclosed, will eventually help illuminate and interpret the hydro-carbonate targets for the seismic processing.

A method and apparatus for imaging seismic data includes obtaining an initial model of a subsurface formation, wherein the model includes a plurality of nodes that form at least part of a grid; an initial dip value for the nodes; and a set of origin coordinates for each of the nodes; performing bottom-up ray tracing for each node in the model, resulting in a set of arrival coordinates for each node; identifying a plurality of gathers from the seismic data; for each gather: calculating a set of midpoint coordinates; defining a midpoint vicinity surrounding the set of midpoint coordinates; identifying the nodes having arrival coordinates within the midpoint vicinity; and estimating a unique aperture for each of the gathers based on the respective origin coordinates; storing the estimated apertures in a table; and generating a subsurface volume or image with subsurface reflectors determined with apertures of the respective gathers.

This disclosure describes processes and systems for generating a high-resolution velocity model of a subterranean formation from recorded seismic data gathers obtained in a marine seismic survey of the subterranean formation. A velocity model is computed by iterative FWI using reflections, resolving the velocity field of deep subterranean targets without requiring ultralong offsets. The processes and systems use of an impedance sensitivity kernel to characterize reflections in a modeled wavefield, and then use the reflections to compute a velocity sensitivity kernel that is used to produce low-wavenumber updates to the velocity model. The iterative process is applied in a cascade such that position of reflectors and background velocity are simultaneously updated. Once the low-wavenumber components of the velocity model are updated, the velocity model is used as an input of conventional FWI to introduce missing velocity components (i.e., high-wavenumber) to increase the resolution of the velocity model.

A tangible, non-transitory computer-readable medium configured to store instructions executable by a processor of an electronic device to access a beam migration image of a subsurface target. In addition, the computer-readable medium is configured to store instructions executable by a processor of an electronic device to determine a decomposition criteria based on at least one of subsurface dip inclinations, subsurface dip azimuths, or a combination thereof. Further, the computer-readable medium is configured to store instructions executable by a processor of an electronic device to decompose the beam migration image into a plurality of partial images according to the decomposition criteria to provide various views of the subsurface target. The plurality of partial images are usable by seismic interpreters in exploration for hydrocarbons within the subsurface target.

A method and a system for implementing the method are disclosed wherein the seismic input data and land acquisition input data may be obtained from a non-flat surface, sometimes mild or foothill topography as well as the shot and receiver lines might not necessarily be straight, and often curve to avoid obstacles on the land surface. In particular, the method and system disclosed, decomposes the cross-spread data into sparse common spread beams, then maps those sparse beams into common-spread depth domain, in order to finally stack them to construct the subsurface depth images. The common spread beam migration and processing have higher signal to noise ratio, as well as faster turn-around processing time, for the cross-spread land acquisition over the common-shot or common offset beam migration/processing. The common spread beam migration method and system disclosed, will eventually help illuminate and interpret the hydro-carbonate targets for the seismic processing.

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.