Abstract: Examples of methods and systems are disclosed. The methods include obtaining seismic data regarding a subsurface region of interest, wherein the seismic data comprises time-space pressure data. The methods also include determining, using a seismic processor, a pressure derivative with respect to depth. The methods further include determining, using the seismic processor, a phase-shifted pressure derivative. The methods still further include determining, using the seismic processor, a transformed phase-shifted pressure derivative. The methods also include determining, using the seismic processor, transformed phase-shifted pressure data based, at least in part, on the transformed phase-shifted pressure derivative. The methods further include determining, using the seismic processor, time-space filtered pressure data. The methods still further include determining, using the seismic processor, a first direction wavefield and a second direction wavefield.

Abstract: Methods and systems are disclosed. The method includes determining a seismic acquisition grid for a seismic survey, composed of a first receiver grid, a second receiver grid, and a source grid, deploying a first plurality of seismic receivers according to the locations of the first receiver grid, and a second plurality of seismic receivers according to locations of the second receiver grid, and activating a seismic source at location based on the source grid. The method further includes recording a first seismic dataset using the first plurality of receivers and a second seismic dataset using the second plurality of receivers, and generating a velocity model using a velocity analysis method based on the first seismic dataset. The method still further includes forming a seismic image based on the velocity model and the second seismic dataset; and determining a location of a hydrocarbon reservoir based on the seismic image.

Abstract: A method may include determining various migrated azimuthal dip-angle gathers based on a migration function and seismic data from a seismic survey regarding a geological region of interest. The method may further include determining various partial dip-angle images using the migrated azimuthal dip-angle gathers. The method may further include determining various azimuthal bins that include the partial dip-angle images. The method may further include determining various wavefield components using an independent component analysis (ICA) function and the partial dip-angle images among the azimuthal bins. The method may further include determining a geological feature within the geological region of interest using the wavefield components.

Abstract: Systems and methods for automated first arrival picking are disclosed. The method includes obtaining a seismic dataset composed of a plurality of seismic gathers and determining a pilot for each gather, where the pilot includes a position on an ordinate axis for each seismic trace representing a first arrival. The method continues iteratively until a stopping criterion is met by creating a preconditioned gather using the pilot, determining a differential pilot using global path tracing subject to a constraint and incrementing the pilot using the differential pilot to create a total picked first arrival. Once the stopping criterion has been met, the method further includes determining a final picked first arrival based on the total picked first arrival, determining a seismic velocity model from the final picked first arrival using a tomographic inversion and creating a seismic image using the seismic velocity model and the seismic dataset.

Abstract: A method is described that includes obtaining a multicomponent seismic data set for a subterranean region of interest and determining, using a computer processor, a PP stacked time-domain seismic image and a PS stacked time-domain seismic image from the multicomponent seismic data set. The method further includes transforming a recording-time axis of at least one of the PP stacked time-domain seismic image and the PS stacked time-domain seismic image to produce a pair of coarsely-registered PP and PS seismic images and filtering at least one of the pair to produce a pair of spectrally-matched PP and PS seismic images. Further, the method includes dynamically warping at least one of the pair of spectrally-matched PP and PS seismic images to produce a pair of fully-registered PP and PS seismic images.

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

Abstract: A method and system for forming a seismic image of a subterranean region are disclosed. The method includes determining an initial plan for a seismic survey with a value for each member of a set of acquisition parameters and acquiring a first seismic dataset from a first portion of the seismic survey based on the initial plan. The method further includes transmitting the first seismic dataset to a seismic processor, determining a first seismic image from the first seismic dataset by performing expedited seismic processing and determining a first updated plan for the seismic survey based on the first seismic image and acquiring a second seismic dataset from a second portion of the seismic survey based on the first updated plan. The method still further includes transmitting the second seismic dataset to the seismic processor and determining the seismic based on the first seismic dataset and the second seismic dataset.

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

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

Abstract: Techniques for determining a wellbore drilling path includes identifying input seismic data associated with a subterranean zone that includes a wellbore drilling target. The input seismic data includes primary seismic events and multiple seismic events. The input seismic data is processed to remove the multiple seismic events and at least one of the primary seismic events from the input seismic data. An orthogonalization of the processed input seismic data is performed to recover the at least one primary seismic event into a seismic image of the subterranean zone that excludes at least a portion of the multiple seismic events. A wellbore path is determined from a terranean surface toward the wellbore drilling target for a drilling geo-steering system based on the seismic image of the subterranean zone.

Abstract: Methods and systems for determining an interbed multiple attenuated pre-stack seismic dataset are disclosed. The methods include forming a post-stack seismic image composed of post-stack traces from the pre-stack seismic dataset and identifying a first, second, and third post-stack horizon on each of the post-stack traces. The methods further include for each pre-stack trace, generating a first, second, and third multiple-generator trace based on the first, second and third post-stack horizon and determining a correlation trace based, at least in part, on a correlation between the first multiple-generator trace and the second multiple-generator trace.

Abstract: Methods and systems, including computer programs encoded on a computer storage medium can be used for identifying primary-wave (P-wave) and secondary-wave (S-wave) characteristics of an underground formation by separating P-wave and S-wave modes of seismic data generated by applying a seismic source to a subterranean region of a geological area. Particle motion vectors of a P-wave are parallel to a propagation vector of the P-wave, whereas particle motion vectors of an S-wave are perpendicular to a propagation vector of the S-wave. The parallel and perpendicular relationship between the motion and propagation vectors of the respective P- and S-waves provide a basis for separating P- and S-wave components from a wavefield. The separation methodology extracts P-wave components and S-wave components from the wavefield based on an estimated angle between propagation vectors and wave motion vectors for the wavefield.

Abstract: A system for seismic imaging of a subterranean geological formation, the system includes a receiver configured to obtain seismic data comprising a data volume representing a post-stacked image. The system includes a filtering module configured to: apply frequency-wavenumber (F-K) filter to the data volume extract a negative-dip structure image and apply the frequency-wavenumber (F-K) filter to the data volume extract a positive-dip structure image. The system includes a diffraction rendering module configured to: multiply the positive-dip structure image with the negative-dip structure image and generate a diffraction-enhanced seismic image representing a geological formation of the data volume.

Abstract: A method is described that includes obtaining a multicomponent seismic data set for a subterranean region of interest and determining, using a computer processor, a PP stacked time-domain seismic image and a PS stacked time-domain seismic image from the multicomponent seismic data set. The method further includes transforming a recording-time axis of at least one of the PP stacked time-domain seismic image and the PS stacked time-domain seismic image to produce a pair of coarsely-registered PP and PS seismic images and filtering at least one of the pair to produce a pair of spectrally-matched PP and PS seismic images. Further, the method includes dynamically warping at least one of the pair of spectrally-matched PP and PS seismic images to produce a pair of fully-registered PP and PS seismic images.

Abstract: A computer-implemented method includes the following. A frequency sweep using sweep parameters is emitted from a vibratory seismic source into geological layers. The sweep parameters include frequencies and modulation parameters for seismic waves. Signals are received from one or more sensors. The signals include seismic data acquisition information, including values identifying energy reflected back from boundaries where rock properties change. A determination is made regarding which of the reflected seismic waves are attenuated. The determination uses an integral transform and a thresholding algorithm for image segmentation. Optimum sweep parameters are determined based on the reflected seismic values that are attenuated and updated to compensate for local geology effects. The emitting, receiving, determining attenuation, determining optimum parameters, and updating are repeated until the received signals are determined to be satisfactory.

Abstract: A computer system receives a post-stack time-domain image having a first spectrum and representing one or more subsurface structures. The computer system reconstructs an increased-frequency version of the post-stack time-domain image using L0-constrained inversion and a least-squares mismatch ratio. The increased-frequency version of the post-stack time-domain image includes structural artifacts. The computer system removes the structural artifacts from the increased-frequency version of the post-stack time-domain image using singular value decomposition. The computer system combines the increased-frequency version of the post-stack time-domain image with the post-stack time-domain image using a weighting function. The computer system generates a combined version of the increased-frequency version of the post-stack time-domain image and the post-stack time-domain image. The combined version represents the one or more subsurface structures and has a second spectrum broader than the first spectrum.

Abstract: A system for seismic imaging of a subterranean geological formation uses a two-way imaging condition. A seismic signal is emitted into a subterranean formation and recorded at receiver(s). Source and receiver wavefields are decomposed into respective right-down/left-up and left-down/right-up propagating waves. The right-down/left-up and left-down/right-up direction can be defined along the direction emitted from the source or receiver to corresponding direction in two dimensional (2D) case. An imaging condition for generating both a positive-dip structure image and a negative-dip structure image is the inner product of the wavefields. Applying the sample-by-sample multiplication imaging condition to the opposite dip images, the diffraction energy is retained while the reflection energy is significantly attenuated. The diffraction image can be used to detect faults and fractures in subsurface regions.

Abstract: The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for suppressing noises in seismic data. One computer-implemented method includes receiving, at a data processing apparatus, a set of seismic data associated with a subsurface region; flattening, by the data processing apparatus, the set of seismic data according to an identified seismic event; dividing, by the data processing apparatus, the set of seismic data into a plurality of spatial windows; randomizing, by the data processing apparatus, the set of seismic data according to a random sequential order; filtering, by the data processing apparatus, the randomized seismic data; and reorganizing, by the data processing apparatus, the filtered seismic data according to a pre-randomization order.

Abstract: A method of generating diffraction images based on wave equations includes generating a source wavefield and a receiver wavefield. Based on the source wavefield, a first source wavefield propagating in a first direction and a second source wavefield propagating in a second direction are generated. Based on the receiver wavefield, a first receiver wavefield propagating in the first direction and a second receiver wavefield propagating in the second direction are generated. A first seismic image is generated based on the first source wavefield and the first receiver wavefield. A second seismic image is generated based on the second source wavefield and the second receiver wavefield. A final seismic image is generated based on the first seismic image and the second seismic image.

Abstract: Methods, systems, and computer-readable medium to perform operations for suppressing seafloor geophone seismic data noise. A computing system applies a vertical geophone de-noise process to enhance a compressional wave signal that is free from (independent of) shear energy leakage. This enhances the signal to noise (S/N) ratio of the vertical geophone component and concurrently make the vertical geophone component consistent with a hydrophone component.