Landmark Graphics Corporation

Atlantic Richfield Company

Halliburton Energy Services, Inc.

Landmark Graphics Corporation, a Halliburton Company

Fleming, Graham C.

Graham Christopher Fleming

Terry Wayne Wong

Qinghua Wang

Qin Lu

James William Watts

Richard Edward Hinkley

Dominic Camilleri

Kefei Wang

Amit Kumar

Shahin Abasov

Ron Mossbarger

Mehmet Tarman

Wei David Liu

Jie Ma

James William Iii Watts

Systems, methods, and computer-readable media for forecasting a time to complete a simulation of a reservoir simulation model. A processor can run a simulation of a reservoir simulation model for a first predetermined period of time associated with a first computation step of the simulation. The processor can measure a performance result of the simulation for the first predetermined period of time. The performance result can include an elapsed wall clock time and a timespan simulated. The processor can then forecast a remaining time to complete the simulation of the reservoir simulation model based on the performance result.

Time-to-finish simulation forecaster

A method for determining active constraint equations in a network of wells and surface facilities includes constructing at least one constraint equation for a connection in the network. Each constraint equation includes a respective slack variable and a respective slack variable multiplier. The method further includes constructing a base equation for the connection. The base equation includes the respective slack variable and another respective slack variable multiplier. The method further includes introducing a pseudo slack variable for another connection in the network such that a Schur complement, of a matrix of constraint and base equations dependent only on slack variable multipliers, is sparse. The method further includes solving for each respective slack variable using the Shur complement matrix. The method further includes adjusting a variable parameter of the network using results from solving for each respective slack variable.

Determining active constraints in network using pseudo slack variables

In some embodiments, a system, as well as a method and an article, may operate to generate map values for a plurality of parameters corresponding to respective grid blocks of a reservoir, wherein the values have been previously generated based on an initial simulation result from a model of the reservoir; to generate a sector surrogate model that includes a subset of grid blocks of the reservoir based on a criterion for identifying grid blocks that negatively affect simulation of the reservoir; to provide data inputs to execute a simulation of the reservoir using the sector surrogate model; and to generate revised data inputs, based on results of the simulation of the reservoir using the sector surrogate model, to use in a subsequent simulation using the model of the reservoir. Additional apparatus, systems, and methods are disclosed.

Model tuning using boundary flux sector surrogates

In some embodiments, a method for locally lumped equation of state fluid characterization can include determining a set of components for the material balance calculations for a plurality of grid blocks of a reservoir. The plurality of grid blocks can experience different recovery methods between them. Lumping schemes can be determined for the plurality of grid blocks. Phase behavior calculations can be performed on the plurality of grid blocks, wherein different lumping schemes can be used across the plurality of grid blocks.

Locally lumped equation of state fluid characterization in reservoir simulation

Systems and methods for determining operating settings of a fluid production system with a gathering network and a plurality of reservoirs, at least some of which include a method that includes lumping common pseudo-components into reservoir pseudo-components associated with each of the reservoirs. Each common pseudo-component is associated with the gathering network and maps to one reservoir pseudo-component of each of the reservoirs. The method further includes performing at least part of a fully-coupled simulation of the gathering network and the reservoirs using the common pseudo-components to obtain one or more operating settings that meet one or more constraints of the production system, and presenting to a user the operating settings.

Systems and methods for solving a multi-reservoir system with heterogeneous fluids coupled to a common gathering network

In some embodiments, a system, as well as a method and an article, may operate to generate a first matrix, based on equations that model a reservoir, that includes mass conservation and volume balance information for grid blocks in the reservoir; to generate a second matrix, based on the first matrix, that includes saturation information and pressure information of each grid block; to remove the saturation information from the second matrix to generate a third matrix that includes only pressure information; to solve the third matrix to generate a first pressure solution; to solve the second matrix based on the first pressure solution to generate a first saturation solution and a second pressure solution; and to use the first saturation solution and the second pressure solution to generate a solution of the first matrix. Additional apparatus, systems, and methods are disclosed.

Multi-stage linear solution for implicit reservoir simulation

A method of managing a network of wells and surface facilities includes partitioning the network to obtain a pipe sub-network and two or more well sub-networks, and constructing an associated set of equations that represent a steady-state fluid flow in the sub-network. The method further includes setting boundary conditions for each well sub-network, and determining a steady-state flow solution for each well sub-network. The method further includes establishing boundary conditions for the pipe sub-network, and finding a steady-state flow solution for the pipe sub-network. If the solution does not match the estimated pressure, the method further includes adjusting the estimated pressure, repeating said setting, determining, establishing, finding, and adjusting operations until the calculated and estimated pressures converge, and analyzing flow rates of the steady-state flow solutions to evaluate suitability of a modification to the network.

Managing a network of wells and surface facilities by finding a steady-state flow solution for a pipe sub-network

System and methods of modeling fluids in a simulation of fluid production in a multi-reservoir system with a common surface network are provided. Pressure-volume-temperature (PVT) data is determined for fluids in each of a plurality of reservoirs coupled to the common surface network. A shared equation of state (EOS) characterization representing each of the fluids across the plurality of reservoirs is generated based on the corresponding PVT data. Data representing properties of the fluids in each reservoir is calculated based on the shared EOS characterization of the fluids. When the calculated data is determined not to match the PVT data associated with the fluids in each reservoir, to the shared EOS characterization is adjusted based on a difference between the calculated data and the PVT data.

Shared equation of state characterization of multiple fluids

System and methods of simulating fluid production in a multi-reservoir system with a common surface network are provided. Black oil data is matched with a common equation of state (EOS) model for each of a plurality of reservoirs. The black oil data representing fluids within each reservoir. At least one multi-dimensional black oil table representing a mix of the fluid components to be produced from each of the plurality of reservoirs via the common surface network is generated based on the EOS model that matches the one or more black oil tables for each reservoir. Properties of the fluids in the mix during a simulation of fluid production from the plurality of reservoirs are calculated based on the generated multi-dimensional black oil table for each reservoir.

Modified black oil model for calculating mixing of different fluids in a common surface network

Oil or gas fields and associated facilities are computer simulated by solving linear algebraic equations using a Schwarz method with at least one well node included in the overlap. According to another aspect, equations for the wellbore nodes are processed to eliminate some non-pressure variables, e.g. component compositions. Other features are also provided.

Oil or gas production using computer simulation of oil or gas fields and production facilities

Systems and methods for determining operating settings for a fluid production system that includes wells and a processing facility, the method including identifying at least one swing well and one or more non-swing wells and determining at least one target parameter equation that is a function of at least one processing facility constraint and one or more non-swing well parameters. The method further includes incorporating the target parameter equation into a fully-coupled equation set representing the production system, simulating the production system using the fully-coupled equation set to obtain one or more operating settings that meet the processing facility constraint, and presenting to a user the operating settings.

Systems and methods for optimizing facility limited production and injection in an integrated reservoir and gathering network

In some embodiments, systems, methods, and articles may operate to compute, in parallel, to determine values of unknowns in network equations associated with a network of sub-surface wells and at least one surface facility, for intra-well subdivisions of the network, and then for inter-well subdivisions of the network, wherein the computing is based on default values of the unknowns, or prior determined values of the unknowns. Additional activities may include constructing a distributed Jacobian matrix having portions comprising coefficients of the unknowns distributed among a number of processors, wherein each of the portions is distributed to a particular one of the processors previously assigned to corresponding ones of the subdivisions. The Jacobian matrix may be factored to provide factors and eliminate some of the unknowns. Back-solving is used to determine remaining unsolved ones of the unknowns, using the factors. Additional apparatus, systems, and methods are described.

Parallel network simulation apparatus, methods, and systems

System and methods of simulating fluid production in a multi-reservoir system with a common surface network are provided. Black oil data is matched with an equation of state (EOS) model representing different fluid components of each reservoir in the multi-reservoir system. The black oil data is converted into a two-component black oil model for each reservoir, based on the EOS model. Fluid production in the multi-reservoir system is simulated for at least one simulation point in the common surface network, based in part on the two-component black oil model of each reservoir. When fluids produced at the simulation point are determined to be from different reservoirs, properties of the fluids are calculated based on weaved EOS models of the different reservoirs. Otherwise, properties of the fluids are calculated using the two-component black oil model for the reservoir from which the fluids are produced.

Simulating fluid production in a common surface network using EOS models with black oil models

Systems and methods for 2D domain decomposition during parallel reservoir simulation to balance the active cells in a reservoir model.

Systems and methods for two-dimensional domain decomposition during parallel reservoir simulation

A disclosed method for a hydrocarbon production system includes collecting production system data. The method also includes performing a simulation based on the collected data, a fluid model, and a fully-coupled set of equations. The method also includes expediting convergence of a solution for the simulation by reducing occurrences of non-physical attributes during the simulation. The method also includes storing control parameters determined for the solution for use with the production system.

Growing community of inventors

Houston, TX, United States of America

Graham Christopher Fleming