Abstract: Improved reservoir simulation methods and systems are provided that employ a new velocity model in conjunction with a sequential implicit (SI) formulation or Sequential Fully Implicit (SF) formulation for solving the discrete form of the system of nonlinear partial differential equations. In embodiments, the new velocity model employs a fluid transport equation part based on calculation of phase velocity for a number of fluid phases that involves capillary pressure and a modification coefficient. In embodiments, the modification coefficient can be based on a derivative of capillary pressure with respect to saturation. In another aspect, the new velocity model can employ an estimate of the phase velocity of the water phase vw_est that is based on one or more derivatives of capillary pressure of the water phase as a function of water saturation.

Abstract: Systems, computer-readable media, and methods for performing a reservoir simulation by obtaining reservoir data; translating the reservoir data into grid properties to create a grid; dividing the grid into domains; generating coarse grids corresponding to each domain; processing the domains, where processing a domain includes: calculating pressure for the domain using a coarse grid corresponding to the domain, calculating flux for the domain using a coarse grid corresponding to the domain, and calculating transport of fluids for the domain using a coarse grid corresponding to the domain; and generating a reservoir simulation corresponding to the grid based on processing each domain. The domains can be processed in parallel on different computer systems, different processors, or different cores.

Abstract: 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.

Abstract: Improved reservoir simulation methods and systems are provided that employ a new velocity model in conjunction with a sequential implicit (SI) formulation or Sequential Fully Implicit (SF) formulation for solving the discrete form of the system of nonlinear partial differential equations. In embodiments, the new velocity model employs a fluid transport equation part based on calculation of phase velocity for a number of fluid phases that involves capillary pressure and a modification coefficient. In embodiments, the modification coefficient can be based on a derivative of capillary pressure with respect to saturation. In another aspect, the new velocity model can employ an estimate of the phase velocity of the water phase ?w_est that is based on one or more derivatives of capillary pressure of the water phase as a function of water saturation.

Abstract: Methods, computing systems, and computer-readable media for multi-scale modeling. The method includes determining a first matrix for a plurality of fine cells of a model based at least in part on a physical property value represented by respective fine cells, identifying one or more overlapped cells of the plurality of fine cells that are part of at least two of the plurality of subdomains, and determining a second matrix. Determining the second matrix includes determining an intermediate product by multiplying the first matrix by a prolongation matrix, which includes predicting a row of zeros in the intermediate product for the plurality of fine cells that are not the one or more overlapped cells and are not part of the at least two of the plurality of subdomains that include the one or more overlapped cells. Determining the second matrix also includes multiplying the intermediate product by a restriction matrix.

Abstract: Systems, computer-readable media, and methods for performing a reservoir simulation by obtaining reservoir data; translating the reservoir data into grid properties to create a grid; dividing the grid into domains; generating coarse grids corresponding to each domain; processing the domains, where processing a domain includes: calculating pressure for the domain using a coarse grid corresponding to the domain, calculating flux for the domain using a coarse grid corresponding to the domain, and calculating transport of fluids for the domain using a coarse grid corresponding to the domain; and generating a reservoir simulation corresponding to the grid based on processing each domain. The domains can be processed in parallel on different computer systems, different processors, or different cores.

Abstract: 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.

Abstract: Computing systems, methods, and computer-readable media for modeling behavior of at least one fluid in a reservoir are disclosed. More particularly, the techniques provide consistent and robust numerical formulations for solutions to linear system of equations arising from the linearization of coupled nonlinear hyperbolic/parabolic (elliptic) partial differential equations (PDEs) of fluid flow in heterogeneous anisotropic porous media.

Abstract: Methods, computing systems, and computer-readable media for multi-scale modeling. The method includes determining a first matrix for a plurality of fine cells of a model based at least in part on a physical property value represented by respective fine cells, identifying one or more overlapped cells of the plurality of fine cells that are part of at least two of the plurality of subdomains, and determining a second matrix. Determining the second matrix includes determining an intermediate product by multiplying the first matrix by a prolongation matrix, which includes predicting a row of zeros in the intermediate product for the plurality of fine cells that are not the one or more overlapped cells and are not part of the at least two of the plurality of subdomains that include the one or more overlapped cells. Determining the second matrix also includes multiplying the intermediate product by a restriction matrix.

Abstract: To simulate a subterranean structure having fracture corridors, a model is used to represent the subterranean structure, where the model also provides a representation of the fracture corridors. A streamline simulation is performed using the model.

Abstract: To simulate a subterranean structure having fracture corridors, a model is used to represent the subterranean structure, where the model also provides a representation of the fracture corridors. A streamline simulation is performed using the model.

Abstract: A method of determining fluid flow in a volume containing two or more fluid components, including determining a pressure field for the volume. One or more streamlines are determined from the pressure field, and the fluid composition is solved for the fluid composition along the, or each, streamline. The pressure may also be solved along the, or each, streamline. The step of solving along the, or each, streamline may be performed using a finite difference technique.

Abstract: A method of determining fluid flow in a volume containing two or more fluid components, comprising determining a pressure field for the volume. One or more streamlines are determined from the pressure field, and the fluid composition is solved for the fluid composition along the or each streamline. The pressure may also be solved along the or each streamline. The step of solving along the or each streamline may be performed using a finite difference technique.