Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media

Ahmad Abushaikha, M. J. Blunt, O. R. Gosselin

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

Abstract

Modelling fluid flow in highly heterogeneous and fractured reservoirs is a challenging task. These reservoirs typically have a complex structure with large and sharp variations in their material properties. Node Control Volume Finite Element (NCVFE) has been used to model those types of reservoirs at the fracture scale for the last decade. However, since the control volumes are constructed around the nodes and the material properties are assigned on elements, there is a loss of accuracy and associated fluid smearing when modelling multi-phase flows. We present a new numerical method to improve the modelling of multi-phase fluid flow in these reservoirs, called Interface Control Volume Finite Element (ICVFE). The method drastically decreases the smearing effects observed with other CVFE methods, such as NCVFE, while being mass conservative and numerically consistent. The pressure is computed at the interfaces of elements, and the control volumes are constructed around them, instead of at the element nodes. This assures that a control volume straddles, at most, two elements, which decreases the fluid smearing between neighbouring elements when large variations in their material properties are present. Lowest order Raviart-Thomas vectorial basis functions are used for the pressure calculation, and Lagrange basis functions are used to compute fluxes. The method is a combination of Mixed Hybrid Finite Element (MHFE) and FE methods. Its accuracy and convergence are tested using three dimensional tetrahedral elements to represent heterogeneous and fractured reservoirs. Our new approach is shown to be more accurate than current methods in the literature. Significance • The ICVFE produces less unphysical flows than NCVFE while honouring the material properties of the domain. • It also models more accurate fluid saturation profiles than NCVFE. • The ICVFE method defines the primary variables (pressure and saturation) on the interfaces of elements. Therefore, it computes a high resolution of the primary variables over the finite element mesh (the number of interfaces is larger than the number of elements). This down-scaling is attractive and convenient since the truncated numerical errors decrease with the increase of degrees of freedom, and conventionally this is achieved by refining the mesh.

Original languageEnglish
Title of host publication14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014
PublisherEuropean Association of Geoscientists and Engineers, EAGE
ISBN (Print)9781634391689
Publication statusPublished - 2014
Externally publishedYes
Event14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014 - Catania
Duration: 8 Sep 201411 Sep 2014

Other

Other14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014
CityCatania
Period8/9/1411/9/14

Fingerprint

finite element method
Interfaces (computer)
fluid flow
Porous materials
porous medium
Flow of fluids
Materials properties
Finite element method
Computer simulation
modeling
Fluids
Multiphase flow
Refining
Numerical methods
fluid
Fluxes
saturation
multiphase flow
numerical method
method

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geotechnical Engineering and Engineering Geology
  • Energy Engineering and Power Technology

Cite this

Abushaikha, A., Blunt, M. J., & Gosselin, O. R. (2014). Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media. In 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014 [A13] European Association of Geoscientists and Engineers, EAGE.

Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media. / Abushaikha, Ahmad; Blunt, M. J.; Gosselin, O. R.

14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014. European Association of Geoscientists and Engineers, EAGE, 2014. A13.

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

Abushaikha, A, Blunt, MJ & Gosselin, OR 2014, Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media. in 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014., A13, European Association of Geoscientists and Engineers, EAGE, 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014, Catania, 8/9/14.
Abushaikha A, Blunt MJ, Gosselin OR. Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media. In 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014. European Association of Geoscientists and Engineers, EAGE. 2014. A13
Abushaikha, Ahmad ; Blunt, M. J. ; Gosselin, O. R. / Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media. 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014. European Association of Geoscientists and Engineers, EAGE, 2014.
@inproceedings{7257000b398d436fb4c2e73dd045ee1f,
title = "Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media",
abstract = "Modelling fluid flow in highly heterogeneous and fractured reservoirs is a challenging task. These reservoirs typically have a complex structure with large and sharp variations in their material properties. Node Control Volume Finite Element (NCVFE) has been used to model those types of reservoirs at the fracture scale for the last decade. However, since the control volumes are constructed around the nodes and the material properties are assigned on elements, there is a loss of accuracy and associated fluid smearing when modelling multi-phase flows. We present a new numerical method to improve the modelling of multi-phase fluid flow in these reservoirs, called Interface Control Volume Finite Element (ICVFE). The method drastically decreases the smearing effects observed with other CVFE methods, such as NCVFE, while being mass conservative and numerically consistent. The pressure is computed at the interfaces of elements, and the control volumes are constructed around them, instead of at the element nodes. This assures that a control volume straddles, at most, two elements, which decreases the fluid smearing between neighbouring elements when large variations in their material properties are present. Lowest order Raviart-Thomas vectorial basis functions are used for the pressure calculation, and Lagrange basis functions are used to compute fluxes. The method is a combination of Mixed Hybrid Finite Element (MHFE) and FE methods. Its accuracy and convergence are tested using three dimensional tetrahedral elements to represent heterogeneous and fractured reservoirs. Our new approach is shown to be more accurate than current methods in the literature. Significance • The ICVFE produces less unphysical flows than NCVFE while honouring the material properties of the domain. • It also models more accurate fluid saturation profiles than NCVFE. • The ICVFE method defines the primary variables (pressure and saturation) on the interfaces of elements. Therefore, it computes a high resolution of the primary variables over the finite element mesh (the number of interfaces is larger than the number of elements). This down-scaling is attractive and convenient since the truncated numerical errors decrease with the increase of degrees of freedom, and conventionally this is achieved by refining the mesh.",
author = "Ahmad Abushaikha and Blunt, {M. J.} and Gosselin, {O. R.}",
year = "2014",
language = "English",
isbn = "9781634391689",
booktitle = "14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014",
publisher = "European Association of Geoscientists and Engineers, EAGE",

}

TY - GEN

T1 - Interface Control Volume Finite Element method for modelling fluid flow in heterogeneous porous media

AU - Abushaikha, Ahmad

AU - Blunt, M. J.

AU - Gosselin, O. R.

PY - 2014

Y1 - 2014

N2 - Modelling fluid flow in highly heterogeneous and fractured reservoirs is a challenging task. These reservoirs typically have a complex structure with large and sharp variations in their material properties. Node Control Volume Finite Element (NCVFE) has been used to model those types of reservoirs at the fracture scale for the last decade. However, since the control volumes are constructed around the nodes and the material properties are assigned on elements, there is a loss of accuracy and associated fluid smearing when modelling multi-phase flows. We present a new numerical method to improve the modelling of multi-phase fluid flow in these reservoirs, called Interface Control Volume Finite Element (ICVFE). The method drastically decreases the smearing effects observed with other CVFE methods, such as NCVFE, while being mass conservative and numerically consistent. The pressure is computed at the interfaces of elements, and the control volumes are constructed around them, instead of at the element nodes. This assures that a control volume straddles, at most, two elements, which decreases the fluid smearing between neighbouring elements when large variations in their material properties are present. Lowest order Raviart-Thomas vectorial basis functions are used for the pressure calculation, and Lagrange basis functions are used to compute fluxes. The method is a combination of Mixed Hybrid Finite Element (MHFE) and FE methods. Its accuracy and convergence are tested using three dimensional tetrahedral elements to represent heterogeneous and fractured reservoirs. Our new approach is shown to be more accurate than current methods in the literature. Significance • The ICVFE produces less unphysical flows than NCVFE while honouring the material properties of the domain. • It also models more accurate fluid saturation profiles than NCVFE. • The ICVFE method defines the primary variables (pressure and saturation) on the interfaces of elements. Therefore, it computes a high resolution of the primary variables over the finite element mesh (the number of interfaces is larger than the number of elements). This down-scaling is attractive and convenient since the truncated numerical errors decrease with the increase of degrees of freedom, and conventionally this is achieved by refining the mesh.

AB - Modelling fluid flow in highly heterogeneous and fractured reservoirs is a challenging task. These reservoirs typically have a complex structure with large and sharp variations in their material properties. Node Control Volume Finite Element (NCVFE) has been used to model those types of reservoirs at the fracture scale for the last decade. However, since the control volumes are constructed around the nodes and the material properties are assigned on elements, there is a loss of accuracy and associated fluid smearing when modelling multi-phase flows. We present a new numerical method to improve the modelling of multi-phase fluid flow in these reservoirs, called Interface Control Volume Finite Element (ICVFE). The method drastically decreases the smearing effects observed with other CVFE methods, such as NCVFE, while being mass conservative and numerically consistent. The pressure is computed at the interfaces of elements, and the control volumes are constructed around them, instead of at the element nodes. This assures that a control volume straddles, at most, two elements, which decreases the fluid smearing between neighbouring elements when large variations in their material properties are present. Lowest order Raviart-Thomas vectorial basis functions are used for the pressure calculation, and Lagrange basis functions are used to compute fluxes. The method is a combination of Mixed Hybrid Finite Element (MHFE) and FE methods. Its accuracy and convergence are tested using three dimensional tetrahedral elements to represent heterogeneous and fractured reservoirs. Our new approach is shown to be more accurate than current methods in the literature. Significance • The ICVFE produces less unphysical flows than NCVFE while honouring the material properties of the domain. • It also models more accurate fluid saturation profiles than NCVFE. • The ICVFE method defines the primary variables (pressure and saturation) on the interfaces of elements. Therefore, it computes a high resolution of the primary variables over the finite element mesh (the number of interfaces is larger than the number of elements). This down-scaling is attractive and convenient since the truncated numerical errors decrease with the increase of degrees of freedom, and conventionally this is achieved by refining the mesh.

UR - http://www.scopus.com/inward/record.url?scp=84908193774&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84908193774&partnerID=8YFLogxK

M3 - Conference contribution

SN - 9781634391689

BT - 14th European Conference on the Mathematics of Oil Recovery 2014, ECMOR 2014

PB - European Association of Geoscientists and Engineers, EAGE

ER -