Realistic simulation of structurally complex reservoirs (SCR) is challenging in at least three ways: (1) geological structures must be represented and discretized accurately on vastly different length scales; (2) extreme ranges and discontinuous variations of material properties have to be associated with the discretized structures and accounted for in the computations; and (3) episodic, highly transient and often localized events such as well shut-in have to be resolved adequately within the overall production history, necessitating a highly adaptive resolution of time. To facilitate numerical experiments that elucidate the emergent properties, typical states and state transitions of SCRs, an application programmer interface (API) called complex systems modelling platform (CSMP ++) has been engineered in ANSI/ISO C ++. It implements a geometry and process-based SCR decomposition in space and time, and uses an algebraic multigrid solver (SAMG) for the spatio-temporal integration of the governing partial differential equations. This paper describes a new SCR simulation workflow including a two-phase fluid flow model that is compared with ECLIPSE in a single-fracture flow simulation. Geologically realistic application examples are presented for incompressible 2-phase flow, compressible 3-phase flow, and pressure-diffusion in a sector-scale model of a structurally complex reservoir.

The permeability of the Earth's crust commonly varies over many orders of magnitude. Flow velocity can range over several orders of magnitude in structures of interest that vary in scale from centimeters to kilometers. To accurately and efficiently model multiphase flow in geologic media, we introduce a fully conservative node-centered finite volume method coupled with a Galerkin finite element method on an unstructured triangular grid with a complementary finite volume subgrid. The effectiveness of this approach is demonstrated by comparison with traditional solution methods and by multiphase flow simulations for heterogeneous permeability fields including complex geometries that produce transport parameters and lengths scales varying over four orders of magnitude.