A method may include obtaining various well logs or various core samples regarding a geological region of interest. The method may further include determining various permeability values, various porosity values, and various dolomite volume fraction values regarding the geological region of interest using the well logs or the core samples. The dolomite volume fraction values may correspond to a percentage of dolomite in a total mineral volume. The method may further include determining, using the porosity values, various permeability thresholds corresponding to various predetermined reservoir qualities. The method may further include generating, using the permeability thresholds, the permeability values, and the dolomite volume fraction values, a reservoir model including various dolomite boundaries defining the predetermined reservoir qualities. The method may further include determining a hydrocarbon trap prediction using the reservoir model.

A method of operating a reservoir simulator can include performing a time step of a reservoir simulation using a spatial reservoir model that represents a subterranean environment that includes a reservoir to generate simulation results for a first time where the simulation results include a front defined by at least in part by a gradient at a position between portions of the spatial reservoir model; predicting a position of the front for a subsequent time step for a corresponding second time using a trained machine model; discretizing the spatial reservoir model locally at the predicted position of the front to generate a locally discretized version of the spatial reservoir model; and performing a time step of the reservoir simulation using the locally discretized version of the spatial reservoir model to generate simulation results for the second time.

A method includes receiving well log data for a plurality of wells. A flag is generated based at least partially on the well log data. The wells are sorted into groups based at least partially on the well log data, the flag, or both. A model is built for each of the wells based at least partially on the well log data, the flag, and the groups.

A method includes, in a computer, generating a discretized model of the subsurface formation in space and time. The discretized model comprises at least one physical parameter of the formation and a relationship between the physical parameter and the physical property. For each spatial location and at each time in the discretized model, a time independent solution to the relationship is calculated. A context is defined of a selected number of grid cells surrounding each spatial location. Dimensionality reduction is performed on each context. Each dimensionality reduced context is input into the computer as a separate earth model to train a machine learning system to determine a relationship between the dimensionality reduced context and the physical property. The trained machine learning system is used to estimate the physical property at each spatial location and at each time.

Systems and methods include a method for generating a real-time reservoir porosity log. Historical gas-porosity data is received from previously-drilled and logged wells. The historical gas-porosity data identifies relationships between gas measurements obtained during drilling and reservoir porosity determined after drilling. A gas-porosity model is trained using machine learning and the historical gas-porosity data. Real-time gas measurements are obtained during drilling of a new well. A real-time reservoir porosity log is generated for the new well using the gas-porosity model and real-time gas measurements.

Method for characterizing subterranean formation is described. One method involves simulating a poroelastic pressure response of known fracture geometry utilizing a geomechanical model to generate a simulated poroelastic pressure response. Compiling a database of simulated poroelastic pressure responses. Measuring a poroelastic pressure response of the subterranean formation during a hydraulic fracturing operation to generate a measured poroelastic pressure response. Identifying a closest simulated poroelastic pressure response in the library of simulated poroelastic pressure response. Estimating a geometrical parameter of a fracture or fractures in the subterranean formation based on the closest simulated poroelastic pressure response.

A method may include obtaining fracture image data regarding a geological region of interest. The method may further include determining various fractures in the fracture image data using a first artificial neural network and a pixel-searching process. The method may further include determining a fracture model using the fractures, a second artificial neural network, and borehole image data. The method may further include determining various fracture permeability values using the fracture model and a third artificial neural network. The method may further include determining various matrix permeability values for the geological region of interest using core sample data. The method may further include generating a coarsened grid model for the geological region of interest using a fourth artificial neural network, the matrix permeability values, and the fracture permeability values.

A method for fracture dynamic hydraulic properties estimation and reservoir simulation may include obtaining a first set of images of a first fracture. The method may include obtaining a first set of fracture detections from the first set of images, generating a plurality of numerical calculations based on the first set of fracture detections, and generating a second model based on the plurality of numerical calculations and the first set of fracture detections. The method may further include obtaining a second set of images of a second fracture of a new reservoir, generating a second set of fracture detections of the second fracture, and generating dynamic hydraulic estimations of the second fracture. The method may also include generating a three-dimensional reservoir simulation and determining a plurality of recovery schemes for the new reservoir.

A method and apparatus for generating a fluid saturation model for a subsurface region. One example method generally includes obtaining a model of the subsurface region; for each of a plurality of fluid types: flooding the subsurface region model with the fluid type to generate a flood model; and running a trial petrophysical inversion with the flood model to generate a trial petrophysical model; identifying potential fluid contact regions in the trial petrophysical models; partitioning the subsurface region model at the identified potential fluid contact regions; and constructing the fluid saturation model from the partitioned subsurface region model.

A reservoir simulation platform is provided. The reservoir simulation platform includes a mimetic finite discretization scheme and an operator-based linearization approach. The reservoir simulation system further includes a parallel framework for coupling the mimetic finite discretization scheme and the operator-based linearization approach.