Imperial College London


Faculty of EngineeringDepartment of Earth Science & Engineering

TOTAL Chair in Geological Fluid Mechanics



+44 (0)20 7594 6538m.d.jackson




1.34Royal School of MinesSouth Kensington Campus






BibTex format

author = {Zhang, Z and Geiger, S and Rood, M and Jacquemyn, C and Jackson, M and Hampson, G and De, Carvalho FM and Marques, Machado Silva CC and Machado, Silva JD and Sousa, MC},
doi = {10.2118/182635-PA},
journal = {SPE Journal},
pages = {1946--1962},
title = {A Tracing Algorithm for Flow Diagnostics on Fully Unstructured Grids With Multipoint Flux Approximation},
url = {},
volume = {22},
year = {2017}

RIS format (EndNote, RefMan)

AB - Flow diagnostics is a common way to rank and cluster ensembles of reservoir models depending on their approximate dynamic behavior before beginning full-physics reservoir simulation. Traditionally, they have been performed on corner-point grids inherent to geocellular models. The rapid-reservoir-modeling (RRM) concept aims at fast and intuitive prototyping of geologically realistic reservoir models. In RRM, complex reservoir heterogeneities are modeled as discrete volumes bounded by surfaces that are sketched in real time. The resulting reservoir models are discretized by use of fully unstructured tetrahedral meshes where the grid conforms to the reservoir geometry, hence preserving the original geological structures that have been modeled.This paper presents a computationally efficient work flow for flow diagnostics on fully unstructured grids. The control-volume finite-element method (CVFEM) is used to solve the elliptic pressure equation. The flux field is a multipoint flux approximation (MPFA). A new tracing algorithm is developed on a reduced monotone acyclic graph for the hyperbolic transport equations of time of flight (TOF) and tracer distributions. An optimal reordering technique is used to deal with each control volume locally such that the hyperbolic equations can be computed in an efficient node-by-node manner. This reordering algorithm scales linearly with the number of unknowns.The results of these computations allow us to estimate swept-reservoir volumes, injector/producer pairs, well-allocation factors, flow capacity, storage capacity, and dynamic Lorenz coefficients, which all help approximate the dynamic reservoir behavior. The total central-processing-unit (CPU) time, including grid generation and flow diagnostics, is typically a few seconds for meshes with O (100,000) unknowns. Such fast calculations provide, for the first time, real-time feedback in the dynamic reservoir behavior while models are prototyped.
AU - Zhang,Z
AU - Geiger,S
AU - Rood,M
AU - Jacquemyn,C
AU - Jackson,M
AU - Hampson,G
AU - De,Carvalho FM
AU - Marques,Machado Silva CC
AU - Machado,Silva JD
AU - Sousa,MC
DO - 10.2118/182635-PA
EP - 1962
PY - 2017///
SN - 1930-0220
SP - 1946
TI - A Tracing Algorithm for Flow Diagnostics on Fully Unstructured Grids With Multipoint Flux Approximation
T2 - SPE Journal
UR -
VL - 22
ER -