Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods

M. D. Jackson, J.L.M.A. Gomes, P. Mostaghimi, J.R. Percival, B.S. Tollit, D. Pavlidis, C.C. Pain, A.H. El-Sheikh, A.H. Muggeridge, M.J. Blunt

Research output: Chapter in Book/Report/Conference proceedingConference contribution

19 Citations (Scopus)

Abstract

We present new approaches to reservoir modeling and flow simulation that dispose of the pillar-grid concept that has
persisted since reservoir simulation began. This results in significant improvements to the representation of multi-scale
geological heterogeneity and the prediction of flow through that heterogeneity. The research builds on 20+ years of
development of innovative numerical methods in geophysical fluid mechanics, refined and modified to deal with the unique
challenges associated with reservoir simulation.
Geological heterogeneities, whether structural, stratigraphic, sedimentologic or diagenetic in origin, are represented as
discrete volumes bounded by surfaces, without reference to a pre-defined grid. Petrophysical properties are uniform within
the geologically-defined rock volumes, rather than within grid-cells. The resulting model is discretized for flow simulation
using an unstructured, tetrahedral mesh that honors the architecture of the surfaces. This approach allows heterogeneity over
multiple length-scales to be explicitly captured using fewer cells than conventional corner-point or unstructured grids.
Multiphase flow is simulated using a novel mixed finite element formulation centered on a new family of tetrahedral
element types, PN(DG)-PN+1, which has a discontinuous Nth-order polynomial representation for velocity and a continuous
(order N+1) representation for pressure. This method exactly represents Darcy force balances on unstructured meshes and
thus accurately calculates pressure, velocity and saturation fields throughout the domain. Computational costs are reduced
through (i) automatic mesh adaptivity in time and space and (ii) efficient parallelization. Within each rock volume, the mesh
coarsens and refines to capture key flow processes, whilst preserving the surface-based representation of geological
heterogeneity. Computational effort is thus focused on regions of the model where it is most required.
Having validated the approach against a set of benchmark problems, we demonstrate its capabilities using a number of
test models which capture aspects of geological heterogeneity that are difficult or impossible to simulate conventionally,
without introducing unacceptably large numbers of cells or highly non-orthogonal grids with associated numerical errors.
Our approach preserves key flow features associated with realistic geological features that are typically lost. The approach
may also be used to capture near wellbore flow features such as coning, changes in surface geometry across multiple
stochastic realizations and, in future applications, geomechanical models with fracture propagation, opening and closing.
Original languageEnglish
Title of host publication2013 SPE Reservoir Simulation Symposium
PublisherSPE International
PagesSPE163633
ISBN (Electronic)978-1-61399-233-3
DOIs
Publication statusPublished - 2013
Event2013 SPE Reservoir Simulation Symposium - The Woodlands, Texas, United States
Duration: 18 Feb 201320 Mar 2013

Conference

Conference2013 SPE Reservoir Simulation Symposium
CountryUnited States
CityThe Woodlands, Texas
Period18/02/1320/03/13

Fingerprint

finite element method
modeling
simulation
fluid mechanics
fracture propagation
geological feature
multiphase flow
pillar
model test
rock
numerical method
saturation
geometry
prediction
cost

Cite this

Jackson, M. D., Gomes, J. L. M. A., Mostaghimi, P., Percival, J. R., Tollit, B. S., Pavlidis, D., ... Blunt, M. J. (2013). Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods. In 2013 SPE Reservoir Simulation Symposium (pp. SPE163633). SPE International. https://doi.org/10.2118/163633-MS

Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods. / Jackson, M. D.; Gomes, J.L.M.A.; Mostaghimi, P.; Percival, J.R.; Tollit, B.S.; Pavlidis, D.; Pain, C.C.; El-Sheikh, A.H.; Muggeridge, A.H.; Blunt, M.J.

2013 SPE Reservoir Simulation Symposium. SPE International, 2013. p. SPE163633.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Jackson, MD, Gomes, JLMA, Mostaghimi, P, Percival, JR, Tollit, BS, Pavlidis, D, Pain, CC, El-Sheikh, AH, Muggeridge, AH & Blunt, MJ 2013, Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods. in 2013 SPE Reservoir Simulation Symposium. SPE International, pp. SPE163633, 2013 SPE Reservoir Simulation Symposium, The Woodlands, Texas, United States, 18/02/13. https://doi.org/10.2118/163633-MS
Jackson MD, Gomes JLMA, Mostaghimi P, Percival JR, Tollit BS, Pavlidis D et al. Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods. In 2013 SPE Reservoir Simulation Symposium. SPE International. 2013. p. SPE163633 https://doi.org/10.2118/163633-MS
Jackson, M. D. ; Gomes, J.L.M.A. ; Mostaghimi, P. ; Percival, J.R. ; Tollit, B.S. ; Pavlidis, D. ; Pain, C.C. ; El-Sheikh, A.H. ; Muggeridge, A.H. ; Blunt, M.J. / Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods. 2013 SPE Reservoir Simulation Symposium. SPE International, 2013. pp. SPE163633
@inproceedings{c9c4ab0cd98744f0a1f8f50061c7c237,
title = "Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods",
abstract = "We present new approaches to reservoir modeling and flow simulation that dispose of the pillar-grid concept that has persisted since reservoir simulation began. This results in significant improvements to the representation of multi-scale geological heterogeneity and the prediction of flow through that heterogeneity. The research builds on 20+ years of development of innovative numerical methods in geophysical fluid mechanics, refined and modified to deal with the unique challenges associated with reservoir simulation. Geological heterogeneities, whether structural, stratigraphic, sedimentologic or diagenetic in origin, are represented as discrete volumes bounded by surfaces, without reference to a pre-defined grid. Petrophysical properties are uniform within the geologically-defined rock volumes, rather than within grid-cells. The resulting model is discretized for flow simulation using an unstructured, tetrahedral mesh that honors the architecture of the surfaces. This approach allows heterogeneity over multiple length-scales to be explicitly captured using fewer cells than conventional corner-point or unstructured grids. Multiphase flow is simulated using a novel mixed finite element formulation centered on a new family of tetrahedral element types, PN(DG)-PN+1, which has a discontinuous Nth-order polynomial representation for velocity and a continuous (order N+1) representation for pressure. This method exactly represents Darcy force balances on unstructured meshes and thus accurately calculates pressure, velocity and saturation fields throughout the domain. Computational costs are reduced through (i) automatic mesh adaptivity in time and space and (ii) efficient parallelization. Within each rock volume, the mesh coarsens and refines to capture key flow processes, whilst preserving the surface-based representation of geological heterogeneity. Computational effort is thus focused on regions of the model where it is most required. Having validated the approach against a set of benchmark problems, we demonstrate its capabilities using a number of test models which capture aspects of geological heterogeneity that are difficult or impossible to simulate conventionally, without introducing unacceptably large numbers of cells or highly non-orthogonal grids with associated numerical errors. Our approach preserves key flow features associated with realistic geological features that are typically lost. The approach may also be used to capture near wellbore flow features such as coning, changes in surface geometry across multiple stochastic realizations and, in future applications, geomechanical models with fracture propagation, opening and closing.",
author = "Jackson, {M. D.} and J.L.M.A. Gomes and P. Mostaghimi and J.R. Percival and B.S. Tollit and D. Pavlidis and C.C. Pain and A.H. El-Sheikh and A.H. Muggeridge and M.J. Blunt",
year = "2013",
doi = "10.2118/163633-MS",
language = "English",
pages = "SPE163633",
booktitle = "2013 SPE Reservoir Simulation Symposium",
publisher = "SPE International",

}

TY - GEN

T1 - Reservoir Modeling for Flow Simulation Using Surfaces, Adaptive Unstructured Meshes, and Control-Volume-Finite-Element Methods

AU - Jackson, M. D.

AU - Gomes, J.L.M.A.

AU - Mostaghimi, P.

AU - Percival, J.R.

AU - Tollit, B.S.

AU - Pavlidis, D.

AU - Pain, C.C.

AU - El-Sheikh, A.H.

AU - Muggeridge, A.H.

AU - Blunt, M.J.

PY - 2013

Y1 - 2013

N2 - We present new approaches to reservoir modeling and flow simulation that dispose of the pillar-grid concept that has persisted since reservoir simulation began. This results in significant improvements to the representation of multi-scale geological heterogeneity and the prediction of flow through that heterogeneity. The research builds on 20+ years of development of innovative numerical methods in geophysical fluid mechanics, refined and modified to deal with the unique challenges associated with reservoir simulation. Geological heterogeneities, whether structural, stratigraphic, sedimentologic or diagenetic in origin, are represented as discrete volumes bounded by surfaces, without reference to a pre-defined grid. Petrophysical properties are uniform within the geologically-defined rock volumes, rather than within grid-cells. The resulting model is discretized for flow simulation using an unstructured, tetrahedral mesh that honors the architecture of the surfaces. This approach allows heterogeneity over multiple length-scales to be explicitly captured using fewer cells than conventional corner-point or unstructured grids. Multiphase flow is simulated using a novel mixed finite element formulation centered on a new family of tetrahedral element types, PN(DG)-PN+1, which has a discontinuous Nth-order polynomial representation for velocity and a continuous (order N+1) representation for pressure. This method exactly represents Darcy force balances on unstructured meshes and thus accurately calculates pressure, velocity and saturation fields throughout the domain. Computational costs are reduced through (i) automatic mesh adaptivity in time and space and (ii) efficient parallelization. Within each rock volume, the mesh coarsens and refines to capture key flow processes, whilst preserving the surface-based representation of geological heterogeneity. Computational effort is thus focused on regions of the model where it is most required. Having validated the approach against a set of benchmark problems, we demonstrate its capabilities using a number of test models which capture aspects of geological heterogeneity that are difficult or impossible to simulate conventionally, without introducing unacceptably large numbers of cells or highly non-orthogonal grids with associated numerical errors. Our approach preserves key flow features associated with realistic geological features that are typically lost. The approach may also be used to capture near wellbore flow features such as coning, changes in surface geometry across multiple stochastic realizations and, in future applications, geomechanical models with fracture propagation, opening and closing.

AB - We present new approaches to reservoir modeling and flow simulation that dispose of the pillar-grid concept that has persisted since reservoir simulation began. This results in significant improvements to the representation of multi-scale geological heterogeneity and the prediction of flow through that heterogeneity. The research builds on 20+ years of development of innovative numerical methods in geophysical fluid mechanics, refined and modified to deal with the unique challenges associated with reservoir simulation. Geological heterogeneities, whether structural, stratigraphic, sedimentologic or diagenetic in origin, are represented as discrete volumes bounded by surfaces, without reference to a pre-defined grid. Petrophysical properties are uniform within the geologically-defined rock volumes, rather than within grid-cells. The resulting model is discretized for flow simulation using an unstructured, tetrahedral mesh that honors the architecture of the surfaces. This approach allows heterogeneity over multiple length-scales to be explicitly captured using fewer cells than conventional corner-point or unstructured grids. Multiphase flow is simulated using a novel mixed finite element formulation centered on a new family of tetrahedral element types, PN(DG)-PN+1, which has a discontinuous Nth-order polynomial representation for velocity and a continuous (order N+1) representation for pressure. This method exactly represents Darcy force balances on unstructured meshes and thus accurately calculates pressure, velocity and saturation fields throughout the domain. Computational costs are reduced through (i) automatic mesh adaptivity in time and space and (ii) efficient parallelization. Within each rock volume, the mesh coarsens and refines to capture key flow processes, whilst preserving the surface-based representation of geological heterogeneity. Computational effort is thus focused on regions of the model where it is most required. Having validated the approach against a set of benchmark problems, we demonstrate its capabilities using a number of test models which capture aspects of geological heterogeneity that are difficult or impossible to simulate conventionally, without introducing unacceptably large numbers of cells or highly non-orthogonal grids with associated numerical errors. Our approach preserves key flow features associated with realistic geological features that are typically lost. The approach may also be used to capture near wellbore flow features such as coning, changes in surface geometry across multiple stochastic realizations and, in future applications, geomechanical models with fracture propagation, opening and closing.

U2 - 10.2118/163633-MS

DO - 10.2118/163633-MS

M3 - Conference contribution

SP - SPE163633

BT - 2013 SPE Reservoir Simulation Symposium

PB - SPE International

ER -