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.

A method for determining an internal multiple attenuated seismic image is disclosed. The method includes obtaining a seismic dataset composed of a plurality of seismic traces and for each seismic trace determining an internal multiple trace based, at least in part, on a nested truncated correlation and a bounded convolution of the seismic trace with itself. The method further includes determining an internal multiple attenuated seismic trace based, at least in part, on subtracting the internal multiple trace from the seismic trace and combining the internal multiple attenuated seismic trace to form the internal multiple attenuated seismic image. A system including a seismic source, a plurality of seismic receivers, and a seismic processor for executing the method is disclosed.

A sensor streamer stretch section. At least some of the example embodiments are methods including measuring at least one parameter related to noise while towing the sensor streamer through a body of water with a towing vessel, and adjusting at least one of a spring constant and a damping coefficient of a stretch section disposed proximate the sensor streamer such that the measured parameter is minimized.

The present invention relates to a one-way wave equation prestack depth migration method using an elastic generalized-screen (EGS) wave propagator capable of efficiently expressing the movement of an elastic wave passing through a mutual mode conversion between a P-wave and an S-wave while propagating boundary surfaces of an underground medium, by expanding, to an elastic wave equation, a conventional scalar generalized-screen (SGS) technique capable of quickly calculating the propagation of a wave in a medium in which there is a horizontal speed change, and according to the present invention, provided is a prestack EGS migration method for seismic wave multi-component data, which: can calculate a wave field with higher accuracy in a medium having a complex structure by expanding up to a second term of a Taylor series expansion of a vertical slowness term of a propagator; includes a mode separation operator in the propagator so as to directly use a shot gather as a migration input, without the need to separate multi-component data into a P-wave and an S-wave, enabling P-wave and S-wave image sections to be generated; and is configured to improve the quality of an S-wave migration image by correcting a polarity conversion in a wave number-frequency domain prior to S-wave imaging.

A method for estimating surface waves generates incident, back-scattered, virtual back-scattered and super-virtual back-scattered traces. The stacked super-virtual back-scattered traces are an estimate of the surface waves.

A marine streamer spread for acquiring seismic data, the spread including a streamer having a first portion and a second portion; the first portion including both first pressure sensors and pressure derivative sensors for acquiring the seismic data; and the second portion including second pressure sensors. The first portion imparts ghost diversity to the seismic data by having a variable-depth profile and the pressure derivative sensors impart polarity diversity to the seismic data.

Disclosed is a system and method for predicting internal multiples generators to correct near surface statics, by estimating a first timing or position associated with reflectors using internal multiple generators identified based on predictive deconvolution operators, estimating a second timing or position associated with the reflectors using the acquired seismic surface data, comparing the first timing or position with the second timing or position for each of the reflectors to determine a travel time delay associated with the reflectors, and correcting the acquired seismic surface data using the travel time delay.

Apparatus, computer instructions and method for deghosting seismic data related to a subsurface of a body of water. The method includes inputting data recorded by detectors that are towed by a vessel, the data being associated with waves travelling from the subsurface to the detectors; applying a migration procedure to the data to determine a first image of the subsurface; applying a mirror migration procedure to the data to determine a second image of the subsurface; joint deconvoluting the first image and the second image for deghosting a reflectivity of the subsurface; and generating a final image of the subsurface based on the deghosted reflectivity of the joint deconvoluting step.

The present disclosure includes a method including determining a spatial region for analysis and selecting a segment of time for analysis, analyzing and correcting a plurality of traces from a plurality of receivers using an iterative non-linear inversion algorithm, wherein each iteration of the non-linear algorithm corrects the plurality of traces using at least one set of parameters defining a microseismic event, determining whether a final stack value of the plurality of traces corrected based on the at least one set of parameters of a final iteration of the iterative non-linear inversion algorithm exceeds a predetermined threshold and upon a determination that the final stack value exceeds the predetermined threshold, detecting a microseismic event defined by the at least one set of parameters of final iteration. The present disclosure also includes associated systems and computer-readable media.