Abstract: Systems and methods for estimating reservoir productivity as a function of position in a subsurface volume of interest are disclosed. Exemplary implementations may: obtain an initial depth uncertainty model; obtain training depth uncertainty parameter values from the non-transient storage medium; obtain corresponding training depth uncertainty values; generate a trained depth uncertainty model by training the initial depth uncertainty model using the training depth uncertainty parameter values and the corresponding training depth uncertainty values; and store the trained depth uncertainty model.

Abstract: Systems and methods are disclosed for associating a standard numerical age to an attribute of geological data from disparate locations. Exemplary implementations may include generating a standardized geological age dataset by standardizing geological data to a global reference age based on a dimension of the geological data, a local geotemporal marker, and a dimension to age function; and storing the standardized geological age dataset.

Abstract: Embodiments of generating a field development plan for a hydrocarbon field development are provided herein. One embodiment comprises generating a plurality of training reservoir models of varying values of input channels of a reservoir template; normalizing the varying values of the input channels to generate normalized values of the input channels; constructing a policy neural network and a value neural network that project a state represented by the normalized values of the input channels to a field development action and a value of the state respectively; and training the policy neural network and the value neural network using deep reinforcement learning on the plurality of training reservoir models with a reservoir simulator as an environment such that the policy neural network generates a field development plan. A field development plan may be generated for a target reservoir on the reservoir template using the trained policy network and the reservoir simulator.

Abstract: Systems and methods for estimating reservoir productivity as a function of position in a subsurface volume of interest are disclosed. The systems and methods may: obtain subsurface data and well data corresponding to a subsurface volume of interest; obtain a parameter model; use the subsurface data and the well data to generate multiple production parameter maps; apply the parameter model to the multiple production parameter maps to generate refined production parameter values; generate multiple refined production parameter graphs; display the multiple refined production parameter graphs; generate one or more user input options; receive a defined well design and the one or more user input options selected by a user to generate limited production parameter values; generate a representation of estimated reservoir productivity as a function of position in the subsurface volume of interest using the defined well design and visual effects; and display the representation.

Abstract: Systems and methods are disclosed for associating a standard numerical age to an attribute of geological data from disparate locations. Exemplary implementations may include generating a standardized geological age dataset by standardizing geological data to a global reference age based on a dimension of the geological data, a local geotemporal marker, and a dimension to age function; and storing the standardized geological age dataset.

Abstract: A method is described for predicting permeability including receiving a 3-D earth model including a volume of interest; generating 2-D property images; receiving 2-D fracture images; training a physics-guided neural network using the 2-D fracture images; and predicting permeability using the physics-guided neural network applied to the 2-D property images. The method is executed by a computer system.

Abstract: Systems, methods, and storage media for detecting a given subsurface event in a subsurface volume of interest are disclosed. Exemplary implementations may: receive a subsurface fiber optic data set; receive a sensor data set using one or more sensors; constrain the subsurface fiber optic data set based on a given parameter value of a given parameter within a certain range to generate a constrained subsurface fiber optic data set; use sets of models to refine the constrained subsurface fiber optic data set to generate a refined subsurface fiber optic data set; estimate an event location of the given subsurface event based on the refined subsurface fiber optic data set; and estimate an origin time based on the event location.

Abstract: A method is described for improving subsalt seismic imaging using wavefield redatuming. The method redatums the source wavefield and receiver wavefield to the salt surface. It cross-correlates the redatumed wavefields and uses the time cross-correlations to update the salt surface. The method may be executed by a computer system.

Abstract: A method is described for high-resolution seismic imaging of complex subsurface volumes using a two-stage least-squares reverse time migration with two objective functions. A traveltime misfit objective function enables long-wavelength imaging and an amplitude misfit objective function enables short-wavelength imaging. The method may be executed by a computer system.

Abstract: A system and method for performing a target-oriented reverse time migration may include obtaining a seismic dataset and a geologic model representative of the subsurface volume of interest; identifying a target based on the seismic dataset and the geologic model; calculating a normal to represent a structural dip direction for at least one point on the target; upward propagating a cluster of beams from the at least one point on the target to a surface representing an acquisition surface of the seismic data; identifying valid source-receiver beam pairs based on coverage of the seismic dataset; accumulating coverage of the valid source-receiver beam pairs; identifying traces from the seismic dataset that fall within the accumulated beam pair coverage; and performing a target-oriented reverse time migration using the identified traces to produce a seismic image.

Abstract: A method is described for high-resolution seismic imaging of complex subsurface volumes using a two-stage least-squares reverse time migration with two objective functions. A traveltime misfit objective function enables long-wavelength imaging and an amplitude misfit objective function enables short-wavelength imaging. The method may be executed by a computer system.

Abstract: A method is described for improving subsalt seismic imaging using wavefield redatuming. The method redatums the source wavefield and receiver wavefield to the salt surface. It cross-correlates the redatumed wavefields and uses the time cross-correlations to update the salt surface. The method may be executed by a computer system.

Abstract: A method is described for seismic imaging of complex subsurface volumes including an optimized partial stack based on ray coverage determined for vector image partition tiles. The method may be executed by a computer system.

Abstract: A system and method for performing a target-oriented reverse time migration may include obtaining a seismic dataset and a geologic model representative of the subsurface volume of interest; identifying a target based on the seismic dataset and the geologic model; calculating a normal to represent a structural dip direction for at least one point on the target; upward propagating a cluster of beams from the at least one point on the target to a surface representing an acquisition surface of the seismic data; identifying valid source-receiver beam pairs based on coverage of the seismic dataset; accumulating coverage of the valid source-receiver beam pairs; identifying traces from the seismic dataset that fall within the accumulated beam pair coverage; and performing a target-oriented reverse time migration using the identified traces to produce a seismic image.

Abstract: Beamed data can be obtained from shot gather data, and debeamed data can be obtained from beamed data. Shot gather data for a geophysical volume of interest in a first domain and/or beamed data in a third domain may be received. The first domain has a component relating to seismic wave propagation time and a spatial component relating to lateral spacing. The third domain has a component relating to local plane wave arrival time at a beam center surface location and a component relating to plane wave arrival directions. Data may be transformed between the first domain and a second domain. The second domain has a frequency component and a spatial component corresponding to the spatial component of the first domain. A forward or reverse transform may be applied to transform shot gather data in the first domain to/from beamed data in the third domain.

Abstract: A computer-implemented method and a system for computational acceleration of seismic data processing are described. The method includes defining a specific non-uniform memory access (NUMA) scheduling for a plurality of cores in a processor according to data to be processed; and running two or more threads through each of the plurality of cores.

Abstract: Beamed data can be obtained from shot gather data, and debeamed data can be obtained from beamed data. Shot gather data for a geophysical volume of interest in a first domain and/or beamed data in a third domain may be received. The first domain has a component relating to seismic wave propagation time and a spatial component relating to lateral spacing. The third domain has a component relating to local plane wave arrival time at a beam center surface location and a component relating to plane wave arrival directions. Data may be transformed between the first domain and a second domain. The second domain has a frequency component and a spatial component corresponding to the spatial component of the first domain. Data may be weighted and/or deweighted in the second domain using a weighting function. A forward transform may be applied to weighted shot gather data in the first domain to obtain beamed data in the third domain.