A computer-implemented method can include the following. Seismic vibratory waves through the Earth along a selected vector path are received. An initial value is selected for a first arrival for each of the seismic vibratory waves. Initial values are determined for travel times and velocities of the seismic vibratory waves. Reversed signs of amplitudes of the seismic vibratory waves are determined and corrected. Time intervals are determined based on the initial values of the travel times. Time windowing and filtering in a frequency domain are performed. Final values are determined for first arrivals and travel times for each of the seismic vibratory waves based on the time windowing and filtering. Final values are determined for velocities of the seismic vibratory waves.

Processes and systems are described for generating an image of a subterranean formation from seismic data recorded during a marine survey that employed a moving vibrational source. Processes and systems compute an up-going pressure wavefield from pressure data and vertical velocity data recorded in the marine survey. A direct incident downgoing vertical velocity wavefield that includes Doppler effects created by the moving vibrational source and characterizes a source wavefield and source ghost of the moving vibrational source is computed and deconvolved from the upgoing pressure wavefield to generate a subsurface reflectivity wavefield. The subsurface reflectivity wavefield is effectively free of contamination from the source wavefield, the source ghost, and the Doppler related effects. Processes and systems generate an image of the subterranean formation based on the subsurface reflectivity wavefield, thereby enhancing resolution of the image by attenuating the source-motion effects, source signature, and source ghost of the moving vibration source.

A method is described for seismic imaging including generating extended image gathers by extended reverse time migration of a seismic dataset using an earth model; processing the extended image gathers to generate processed image gathers; performing extended modeling based on the processed image gathers to generate a modeled seismic dataset; enhancing the processed image gathers to generate an enhanced image; performing extended modeling based on the enhanced image gathers to generate a modeled enhanced dataset; differencing the modeled enhanced dataset and the modeled seismic dataset to determine a data residual; inverting the data residual to generate a model residual; updating the earth model based on the model residual to create an updated earth model; performing seismic imaging of the seismic dataset using the updated earth model to create an improved seismic image. The method may be executed by a computer system.

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.

Methods, apparatus, and systems for mapping surface and near surface features and processing artifacts from stacked and processed seismic data are disclosed. In some implementations, a computer system receives a three-dimensional (3D) seismic cube including seismic reflectivity data obtained at a geographical location. The computer system generates a vertical analysis window of the 3D seismic cube. The computer system extracts a second 3D seismic cube from the seismic reflectivity data based on the vertical analysis window. The second 3D seismic cube has multiple vertical amplitude traces associated with the seismic reflectivity data. The computer system generates 3D data comprising multiple frequency traces from the multiple vertical amplitude traces using a domain transform. The computer system generates a two-dimensional (2D) map from the 3D data. The 2D map represents geographical features of the geographical location.

Systems and methods for training a model that refines estimated parameter values include computer processors and non-transitory electronic storage that stores subsurface map data sets that correspond to different subsurface volumes of interest, the system configured to obtain training data including unrefined subsurface map data sets specifying estimated parameter values of a first parameter as a function of position within corresponding subsurface volumes of interest, obtain an initial seismic mapping model, generate a conditioned seismic mapping model, and store the conditioned seismic mapping model in the electronic storage.

A method and apparatus for hydrocarbon management including generating an image of a subsurface formation by: obtaining simultaneous-source survey data, an earth model, and a first and a second velocity model of the subsurface formation; generating synthetic survey data with at least one of the earth model, the first velocity model, and the second velocity model. The method and apparatus may include directly migrating the simultaneous-source survey data; migrating the synthetic survey data; and subtracting the migrated synthetic survey data from the migrated simultaneous-source survey data. The method and apparatus may include subtracting the synthetic survey data from the simultaneous-source survey data; and directly migrating the result of the subtraction. The method and apparatus may include generating an artifact-reduced image.

A method of at least partially deghosting recorded seismic s-waves, wherein recorded seismic data is provided, wherein the recorded seismic data has been recorded at a receiver located beneath the Earth's surface, and wherein the recorded seismic data includes s-wave data. The method may include the steps of finding a model of the Earth's crust for use in deghosting the recorded seismic data using the s-wave data, wherein the model includes at least one region and wherein the model includes the Earth's surface and the location of the receiver, using the model to find a deghosting operator that, when applied to the s-wave data, at least partially deghosts the s-wave data, and applying the deghosting operator to the s-wave data to at least partially deghost the s-wave data.

Prestack images from the object store are hierarchically combined to generate a hierarchically stacked image. The hierarchically stacked image is generated by combining stacked images that includes a stacked image. The stacked image is generated by combining at least the prestack images. Based at least on the hierarchically stacked image, a quality measure of a prestack image is generated. Prior to deleting at least a subset of the prestack images from the object store and based at least on the quality measure, the prestack images are further combined to generate an enhanced stacked image. The stacked image is substituted using the enhanced stacked image. Subsequent to the substituting and prior to deleting at least the subset of the stacked images from the object store, the stacked images are combined to generate an enhanced hierarchically stacked image. The enhanced stacked image and the enhanced hierarchically stacked image are generated using failure recovery metadata. The enhanced hierarchically stacked image is presented.

The present disclosure provides a method of high-resolution amplitude-preserving seismic imaging for a subsurface reflectivity model, including: performing reverse time migration (RTM) to obtain an initial imaging result, performing Born forward modeling on the initial imaging result to obtain seismic simulation data, and performing RTM on the seismic simulation data to obtain a second imaging result; performing curvelet transformation on the two imaging results, performing pointwise estimation in a curvelet domain, and using a Wiener solution that matches two curvelet coefficients as a solution of a matched filter; and applying the estimated matched filter to the initial imaging result to obtain a high-resolution amplitude-preserving seismic imaging result.