Title:

Geothermal Modeling in Complex Geological Systems with the ComPASS Code

Authors:

Simon LOPEZ, Roland MASSON, Feng XING, Laurence BEAUDE, Nabil BIRGLE, Farid SMAI, Konstantin BRENNER, Michel KERN, Gabriel COURRIOUX, Severine CARITG, Yannis LABEAU

Key Words:

reservoir modeling, hydrothermal simulations, multiphase flow, unstructured meshes, fracture network

Conference:

Stanford Geothermal Workshop

Year:

2018

Session:

Modeling

Language:

English

Paper Number:

Lopez

File Size:

1352 KB

View File:

Abstract:

Deep geothermal systems often lie in complex geological settings, with multi-scale geological structures that exert a dominant control on convective processes and the transfer of geothermal heat. Methods based on the implicit description of geometrical objects offer an efficient framework to quickly build structural models of such contexts with the occurrence of discontinuities like faults and fractures. Yet, when it comes to discretizing such models the implicit nature of surfaces make volume meshing a non-trivial task and the results are unstructured polyhedral meshes. Over the last few years, much progress has been made towards the consistent and robust discretization of diffusion processes in porous media. These research efforts resulted in several numerical schemes designed with a sound mathematical basis and able to deal with subsurface spatial heterogeneities (permeability variations, anisotropies…) and general polyhedral meshes. We introduce hereafter the ComPASS platform, an open source initiative that aims at building a geothermal simulation platform relying on one such scheme and recent numerical techniques. The current code is able to handle compositional multiphase Darcy flows, relying on a Coats type formulation, coupled to the conductive and convective transfers of energy. Simulations can be run on unstructured meshes including complex networks of fractures with intersecting fractures. Flow inside the fractures is modelled with a so called hybrid-dimensional model, using a 2D model in the fractures that can have variable apertures and permeability and is coupled with 3D transfers in the matrix. The physics is discretized using a fully implicit time integration combined with the Vertex Approximate Gradient (VAG) finite volume scheme. The fully coupled system is assembled and solved in parallel using the PETSc library and can be run on large computing clusters. An efficient preconditioner is implemented to solve the linear systems at each time step and each Newton type iteration of the simulation. A high level interface to describe the simulations is provided by the Python language, whereas the core routine are written in Fortran and C++. This paper reviews the theoretical foundations and presents a few examples of applications along with the current development perspectives.


ec2-54-225-56-41.compute-1.amazonaws.com, you have accessed 0 records today.

Press the Back button in your browser, or search again.

Copyright 2018, Stanford Geothermal Program: Readers who download papers from this site should honor the copyright of the original authors and may not copy or distribute the work further without the permission of the original publisher.


Attend the nwxt Stanford Geothermal Workshop, click here for details.

Accessed by: ec2-54-225-56-41.compute-1.amazonaws.com (54.225.56.41)
Accessed: Friday 29th of March 2024 04:25:07 AM