A method is disclosed that includes: obtaining a seismic data volume for a subterranean region of interest; transforming, by a computer processor using a non-stationary series analysis, the seismic data volume into a seismic spectral volume where the seismic spectral volume includes a seismic spectrum for each of a plurality of voxels; and determining a seismic attribute volume composed of a seismic attribute for each of the plurality of voxels. The seismic attribute for a voxel of the plurality of voxels is based, at least in part, on an integral of the seismic spectrum for the voxel over a range bounded by a first frequency and a second frequency. The method further includes determining a presence of hydrocarbon in the subterranean region of interest based on the seismic attribute volume. A system for performing the method is also disclosed and described.

A method for estimating a thickness of a deep reservoir may include obtaining seismic data relating to the deep reservoir. The method may include performing spectral decomposition to obtain one or more frequency components from the seismic data. The method may include identifying a number of mono-frequency horizons corresponding to high frequencies in the seismic data, determining whether the deep reservoir is a thin reservoir based on the number of mono-frequency horizons, and estimating the thickness of the deep reservoir when the deep reservoir is determined to be the thin reservoir.

A method of processing a three-dimensional (3D) seismic dataset includes: producing a relative amplitude preserved (RAP) processed seismic volume that includes a reservoir from the 3D seismic dataset; decomposing the RAP processed seismic volume into a first iso-frequency volume; generating a first overburden sediment interval map and a first reservoir interval zone map for the first iso-frequency volume; generating and outputting a first ratio map for the first iso-frequency volume by dividing the first overburden sediment interval map by the first reservoir interval zone map; and characterizing a fluid content of the reservoir based on the first ratio map.

A method for training a backpropagation-enabled segmentation process is used for identifying an occurrence of a sub-surface feature. A multi-dimensional seismic data set with an input dimension of at least two is inputted into a backpropagation-enabled process. A prediction of the occurrence of the subsurface feature has a prediction dimension of at least 1 and is at least 1 dimension less than the input dimension.

A method for training a backpropagation-enabled regression process is used for predicting values of an attribute of subsurface data. A multi-dimensional seismic data set with an input dimension of at least two is inputted into a backpropagation-enabled process. A predicted value of the attribute has a prediction dimension of at least 1 and is at least 1 dimension less than the input dimension.

A method can include receiving a first trained machine model trained via unsupervised learning using unlabeled seismic image data; receiving labeled seismic image data acquired via an interactive interpretation process; and building a second trained machine model, as initialized from the first trained machine model, via supervised learning using the received labels, where the second trained machine model predicts stratigraphy of a geologic region from seismic image data of the geologic region.

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.

A method that includes transforming a relative amplitude preserved 3D seismic volume acquired in the time-domain into a plurality of isofrequency volumes, extracting from the plurality of isofrequency volumes a first isofrequency spectral amplitude volume and a second isofrequency spectral amplitude volume. The method further includes determining an attribute volume from the two isofrequency spectral amplitude volumes, and determining a presence of gas in a subterranean region of interest based on the attribute volume.

The present disclosure describes a computer-implemented method that includes: receiving a seismic dataset of a surveyed subsurface of a reservoir, the seismic dataset comprising observed pressure and production data of the reservoir as well as a set of geological and geo-mechanical parameters representing physical features of the surveyed subsurface; generating multiple realizations of a discrete fracture network (DFN) based on a subset of the set of geological and geo-mechanical parameters; selecting, from the multiple realizations, one or more realizations based on a parameter with a value under a 10% quantile of a full range of likely values; performing a forward simulation for the reservoir based on the selected one or more realizations and the observed pressure and production data; determining that a misfit of the forward simulation is below a threshold based on evaluating an objective function; and producing a model of the reservoir based on the forward simulation.

A method of processing a three-dimensional (3D) seismic dataset includes: producing a relative amplitude preserved (RAP) processed seismic volume that includes a reservoir from the 3D seismic dataset; decomposing the RAP processed seismic volume into a first iso-frequency volume; generating a first overburden sediment interval map and a first reservoir interval zone map for the first iso-frequency volume; generating and outputting a first ratio map for the first iso-frequency volume by dividing the first overburden sediment interval map by the first reservoir interval zone map; and characterizing a fluid content of the reservoir based on the first ratio map.