Abstract: A method for history matching a reservoir model based on actual production data from the reservoir over time generates an ensemble of reservoir models using geological data representing petrophysical properties of a subterranean reservoir. Production data corresponding to a particular time instance is acquired from the subterranean reservoir. Normal score transformation is performed on the ensemble and on the acquired production data to transform respective original distributions into normal distributions. The generated ensemble is updated based on the transformed acquired production data using an ensemble Kalman filter (EnKF). The updated generated ensemble and the transformed acquired production data are transformed back to respective original distributions. Future reservoir behavior is predicted based on the updated ensemble.

Abstract: Systems and methods for modeling petroleum reservoir properties using a gridless reservoir simulation model are provided. Data relating to geological properties of a reservoir formation is analyzed. A tiered hierarchy of geological elements within the reservoir formation is generated at different geological scales, based on the analysis. The geological elements at each of the different geological scales in the tiered hierarchy are categorized. Spatial boundaries between the categorized geological elements are defined for each of the geological scales in the tiered hierarchy. A scalable and updateable gridless model of the reservoir formation is generated, based on the spatial boundaries defined for at least one of the geological scales in the tiered hierarchy.

Abstract: A system is described for estimating well production and injection rates of a subterranean reservoir using machine learning models. The system may include a processor and a non-transitory computer-readable medium comprising instructions that are executable by the processor to cause the processor to perform various operations. The processor may receive a set of static geological data about at least one subterranean reservoir in a subterranean formation. The processor may apply a trained convolutional neural network to the set of static geological data and data on initial states of dynamic reservoir properties to determine dynamic outputs of the subterranean reservoir. The processor may determine well data by extracting the set of static geological data and the dynamic outputs at well trajectories. And, the processor may apply a trained artificial neural network to the well data and subterranean grid information about the subterranean reservoir to generate estimated well production and injection rates.

Abstract: A method for identifying a flow parameter in a wellbore may comprise identifying a state vector at a moment t, performing a flow simulation using a flow model, predicting the state vector and a covariance matrix at the moment t, updating the state vector with an EnKF algorithm, correcting the state vector at the moment t, and updating the flow simulation model. A system for identifying a flow parameter in a wellbore may comprise a distributed acoustic system into a wellbore and an information handling system. The distributed acoustic system may comprise a fiber optic cable and at least one measurement device.

Abstract: A method for creating a seamless scalable geological model may comprise identifying one or more geological scales, establishing a geological tied system, identifying one or more graphical resolution levels for each of the one or more geological scales, constructing the seamless scalable geological model, and producing a post-process model. A system for creating a seamless scalable geological model may comprise an information handling system, which may comprise a random access memory, a graphics module, a main memory, a secondary memory, and one or more processors configured to run a seamless scalable geological model software.

Abstract: A method includes collecting a first set of borehole gravity data at a first time step along a length of a first wellbore and collecting a second set of borehole gravity data at the first time step along a length of a second wellbore. The method also includes interpolating a third set of borehole gravity data at the first time step in an area between the first wellbore and the second wellbore using the first and the second sets of borehole gravity data. Further, the method includes determining a first fluid saturation and a fluid saturation change over time in a reservoir containing the first wellbore and the second wellbore using the first set, the second set, and the third set. Moreover, the method includes controlling wellbore production operations or wellbore injection operations at the first wellbore based on the fluid saturation change.

Abstract: A system is described for estimating well production and injection rates of a subterranean reservoir using machine learning models. The system may include a processor and a non-transitory computer-readable medium comprising instructions that are executable by the processor to cause the processor to perform various operations. The processor may receive a set of static geological data about at least one subterranean reservoir in a subterranean formation. The processor may apply a trained convolutional neural network to the set of static geological data and data on initial states of dynamic reservoir properties to determine dynamic outputs of the subterranean reservoir. The processor may determine well data by extracting the set of static geological data and the dynamic outputs at well trajectories. And, the processor may apply a trained artificial neural network to the well data and subterranean grid information about the subterranean reservoir to generate estimated well production and injection rates.

Abstract: Systems and methods for modeling petroleum reservoir properties using a gridless reservoir simulation model are provided. Data relating to geological properties of a reservoir formation is analyzed. A tiered hierarchy of geological elements within the reservoir formation is generated at different geological scales, based on the analysis. The geological elements at each of the different geological scales in the tiered hierarchy are categorized. Spatial boundaries between the categorized geological elements are defined for each of the geological scales in the tiered hierarchy. A scalable and updateable gridless model of the reservoir formation is generated, based on the spatial boundaries defined for at least one of the geological scales in the tiered hierarchy.

Abstract: A method for creating a seamless scalable geological model may comprise identifying one or more geological scales, establishing a geological tied system, identifying one or more graphical resolution levels for each of the one or more geological scales, constructing the seamless scalable geological model, and producing a post-process model. A system for creating a seamless scalable geological model may comprise an information handling system, which may comprise a random access memory, a graphics module, a main memory, a secondary memory, and one or more processors configured to run a seamless scalable geological model software.

Abstract: A method for history matching a reservoir model based on actual production data from the reservoir over time generates an ensemble of reservoir models using geological data representing petrophysical properties of a subterranean reservoir. Production data corresponding to a particular time instance is acquired from the subterranean reservoir. Normal score transformation is performed on the ensemble and on the acquired production data to transform respective original distributions into normal distributions. The generated ensemble is updated based on the transformed acquired production data using an ensemble Kalman filter (EnKF). The updated generated ensemble and the transformed acquired production data are transformed back to respective original distributions. Future reservoir behavior is predicted based on the updated ensemble.

Abstract: An apparatus for generating forecasts from a high-dimensional parameter data space comprising a reservoir model, a model order reduction module, and an assisted history matching module. The reservoir model having input variables, output variables, and an algorithmic model. The input variables, output variables, and the algorithmic model are generated by a flow simulator module and from a formation and reservoir properties database and a field production database. The model order reduction module generates a subset of the original or transformed input variables. This subset has a reduced parameter space than that of the input variables. The subset is generated using a function decomposition and a design of experiments (sensitivity analysis) to reduce number of original variables and identify original or transformed input variables that can be used to approximate output variables.

Abstract: An apparatus used to generate forecasts from a high-dimensional parameter data space. The apparatus comprising a reservoir model and a flow simulator module. The reservoir model comprising a plurality input variables, output variables, and at least one algorithmic model. The input variables and output variables are generated by the flow simulator module and variables from a formation and reservoir properties database and a field production database. The flow simulator module generates the at least one algorithmic model and the output variables using at least one selected from a group comprising a full-physics flow simulator, proxy flow simulator for assisted history matching, and a proxy flow simulator for field development optimization. The full-physics flow simulator and the two proxy flow simulators generate the at least one algorithmic model using at least one selected from a group comprising the reservoir model, history matching input variables, and optimization input variables.

Abstract: Using production data and a production flow record based on the production data, a deep neural network (DNN) is trained to model a proxy flow simulation of a reservoir. The proxy flow simulation of the reservoir is performed, using an ensemble Kalman filter (EnKF), based on the trained DNN. The EnKF assimilates new data through updating a current ensemble to obtain history matching by minimizing a difference between a predicted production output from the proxy flow simulation and measured production data from a field. Using the updated current ensemble, a second proxy flow simulation of the reservoir is performed based on the trained DNN. The assimilating and the performing are repeated while new data is available for assimilating. Predicted behavior of the reservoir is determined based on the proxy flow simulation of the reservoir. An indication of the predicted behavior is provided to facilitate production of fluids from the reservoir.

Abstract: A method for identifying a flow parameter in a wellbore may comprise identifying a state vector at a moment t, performing a flow simulation using a flow model, predicting the state vector and a covariance matrix at the moment t, updating the state vector with an EnKF algorithm, correcting the state vector at the moment t, and updating the flow simulation model. A system for identifying a flow parameter in a wellbore may comprise a distributed acoustic system into a wellbore and an information handling system. The distributed acoustic system may comprise a fiber optic cable and at least one measurement device.

Abstract: Microseismic-event data can be corrected (e.g., to reduce or eliminate bias). For example, a first distribution of microseismic events that occurred in a first area of a subterranean formation can be determined. The first distribution can be used as a reference distribution. A second distribution of microseismic events that occurred in a second area of the subterranean formation can also be determined. The second area of the subterranean formation can be farther from an observation well than the first area. The second distribution can be corrected by including, in the second distribution, microseismic events that have characteristics tailored for reducing a difference between the second distribution and the first distribution.