A fracture pattern characterization method, system, and computer program product, include classifying fracture patterns of reservoir layers using a correlation of reservoir information associated with target images and the target images extracted from a database and determining a fracture pattern characterization of a new reservoir layer input based on the classified fracture patterns.

The disclosed embodiments include devices and methods to perform inversion processing of well log data. In one embodiment, a method to perform inversion processing of well log data includes obtaining an initial model of an earth formation based on a plurality of modeling parameters that includes formation parameters of the earth formation and calibration factors associated with orientations of antennas of a logging tool utilized to measure raw measurements of the earth formation. The method also includes performing a forward modeling of the modeling parameters to obtain a modeling response, and performing a joint cost function of the first modeling response and raw measurements obtained by the logging tool. The method further includes readjusting the initial model if a result of the joint cost function is not below a threshold, and providing the modeling response if the result of the joint cost function is below the threshold.

The invention pertains to a method for estimating faults in a three-dimensional seismic image block. Directrices are generated within respective first cross-sections of the seismic image block based on points selected by a user. Similarly, generatrices are generated within respective second cross-sections of the block based on points selected by the user. The user inputs relationships between directrices and generatrices. A fault is estimated within the seismic image block as a surface including at least one directric and at least one generatrix having a relationship therebetween.

A method includes receiving historical well-production data, and determining curve parameters based on the historical well production data. Determining the curve parameters includes accounting for uncertainty. The method also includes determining curve fits for the historical data based on the curve parameters using a plurality of models, calculating accuracy for each of the models based on a comparison of the curve fits to the well-production data, comparing the accuracy for each of the models using an information criteria that accounts for uncertainty, selecting one of the models based on the comparison, and determining an estimated ultimate recovery (EUR) for the well using at least the selected one of the models.

A method and system for automating a reservoir simulation. The method includes identifying a simulation parameter associated with a simulation resource to perform a computer-based reservoir simulation using reservoir data associated with a subterranean reservoir and configuring the simulation resource using a simulation engine to include the simulation parameter for performing the reservoir simulation with a reduced likelihood of simulation failure. The method also includes performing the reservoir simulation using the configured simulation resource and the reservoir data to generate reservoir simulation data and evaluate the reservoir.

Embodiments include a method that includes extracting a low resolution formation property model and a high resolution formation property model from at least one log. The method also includes splitting the extracted low resolution formation property model into one or more property compositions. The method also includes generating a high resolution virtual core via the low resolution formation property model and at least a second high resolution formation property, the high resolution virtual core utilizing a moving window analysis to accommodate between different resolutions of the second high resolution formation property model and the low resolution property composition.

A computer-implemented method for updating subsurface models including: using an offset continuation approach to update the model, and at each stage defining a new objective function where a maximum offset for each stage is set, wherein the approach includes, performing a first stage iterative full wavefield inversion with near offset data, as the maximum offset, to obtain velocity and density or impedance models, performing subsequent stages of iterative full wavefield inversion, each generating updated models, relative to a previous stage, wherein the subsequent stages include incrementally expanding the maximum offset until ending at a full offset, wherein a last of the stages yields finally updated models, the subsequent stages use the updated models as starting models, and the full wavefield inversions include constraining scales of the velocity model updates at each stage of inversion as a function of velocity resolution; and using the finally updated models to prospect for hydrocarbons.

The process comprises positioning wells one after another in a group of potential cells of a geocellular model, each positioning of a well comprises: calculating for each cell of the group of potential cells, a fluid property insertion point driver (DFP1) representative of a fluid property maximization; calculating for each cell of the group of potential cells, a maximized distance insertion point driver (DMD1) representative of a maximization of a distance to another cell or group of cells having at least an undesired property; calculating for each cell of the group of potential cells a combined insertion point driver based on the fluid property insertion point driver (DFP1) and the maximized distance insertion point driver (DMD1); defining a well insertion point of the well being positioned at the cell having a maximal combined insertion point driver.

Systems, computer-readable media, and methods are described for performing a reservoir simulation by obtaining reservoir data, obtaining simulation parameters, determining partial differential equations based on the simulation parameters, and performing a timestep of the reservoir simulation based on the reservoir data and the partial differential equations by removing an effect of long coherent structures with high contrasts, such as fractures, faults, high and low permeability channels, or shale layers, from the partial differential equations to generate adapted partial differential equations, and performing an algebraic multiscale method on the adapted partial differential equations to generate an approximated solution. The approximated solution can be used in a subsequent timestep of the reservoir simulation.

Analysis and display of source-to-sink information according to some aspects includes grouping geochronological data associated with a sediment sample into optimized subpopulations within a reference population and target populations, and producing Gaussian functions for the reference population and the target populations using the subpopulations as a priori constraints. The Gaussian functions describe a distribution of zircons. The subpopulations within the reference population and the target populations are compared based on at least one statistical attribute from the Gaussian functions to identify areas of sediment provenance, and the areas of sediment provenance are displayed in various ways, for example, on a paleographic map as of an age of deposition of the sediment sample. A sink-to-sink analysis can also be performed to identify dissimilarities between samples.