-
A Line Search Algorithm for Multiphysics Problems with Fracture Deformation
Authors:
Ivar Stefansson
Abstract:
Models for multiphysics problems often contain strong nonlinearities. Including fracture contact mechanics introduces discontinuities at the transition between open and closed or sliding and sticking fractures. The resulting system of equations is highly challenging to solve. The naïve choice of Newton's method frequently fails to converge, calling for more refined solution techniques such as line…
▽ More
Models for multiphysics problems often contain strong nonlinearities. Including fracture contact mechanics introduces discontinuities at the transition between open and closed or sliding and sticking fractures. The resulting system of equations is highly challenging to solve. The naïve choice of Newton's method frequently fails to converge, calling for more refined solution techniques such as line search methods. When dealing with strong nonlinearities and discontinuities, a global line search based on the magnitude of the residual of all equations is at best costly to evaluate and at worst fails to converge. We therefore suggest a cheap and reliable approach tailored to the discontinuities. Utilising adaptive variable scaling, the algorithm uses a line search to identify the transition between contact states. Then, a solution update weight is chosen to ensure that no fracture cells move too far beyond the transition. We demonstrate the algorithm on a series of test cases for poromechanics and thermoporomechanics in fractured porous media. We consider both single- and multifracture cases and study the importance of proper scaling of variables and equations.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Flexible and rigorous numerical modelling of multiphysics processes in fractured porous media using PorePy
Authors:
Ivar Stefansson,
Jhabriel Varela,
Eirik Keilegavlen,
Inga Berre
Abstract:
Multiphysics processes in fractured porous media is a research field of importance for several subsurface applications and has received considerable attention over the last decade. The dynamics are characterised by strong couplings between processes as well as interaction between the processes and the structure of the fractured medium itself. The rich range of behavior calls for explorative mathem…
▽ More
Multiphysics processes in fractured porous media is a research field of importance for several subsurface applications and has received considerable attention over the last decade. The dynamics are characterised by strong couplings between processes as well as interaction between the processes and the structure of the fractured medium itself. The rich range of behavior calls for explorative mathematical modelling, such as experimentation with constitutive laws and novel coupling concepts between physical processes. Moreover, efficient simulations of the strong couplings between multiphysics processes and geological structures require the development of tailored numerical methods.
We present a modelling framework and its implementation in the open-source simulation toolbox PorePy, which is designed for rapid prototyping of multiphysics processes in fractured porous media. PorePy uses a mixed-dimensional representation of the fracture geometry and generally applies fully implicit couplings between processes. The code design follows the paradigms of modularity and differentiable programming, which together allow for extreme flexibility in experimentation with governing equations with minimal changes to the code base. The code integrity is supported by a multilevel testing framework ensuring the reliability of the code.
We present our modelling framework within a context of thermo-poroelasticity in deformable fractured porous media, illustrating the close relation between the governing equations and the source code. We furthermore discuss the design of the testing framework and present simulations showcasing the extendibility of PorePy, as well as the type of results that can be produced by mixed-dimensional simulation tools.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
Numerical modelling of convection-driven cooling, deformation and fracturing of thermo-poroelastic media
Authors:
Ivar Stefansson,
Eirik Keilegavlen,
Sæunn Halldórsdóttir,
Inga Berre
Abstract:
Convection-driven cooling in porous media influences thermo-poro-mechanical stresses, thereby causing deformation. These processes are strongly influenced by the presence of fractures, which dominate flow and heat transfer. At the same time, the fractures deform and propagate in response to changes in the stress state. Mathematically, the model governing the physics is tightly coupled and must acc…
▽ More
Convection-driven cooling in porous media influences thermo-poro-mechanical stresses, thereby causing deformation. These processes are strongly influenced by the presence of fractures, which dominate flow and heat transfer. At the same time, the fractures deform and propagate in response to changes in the stress state. Mathematically, the model governing the physics is tightly coupled and must account for the strong discontinuities introduced by the fractures. Over the last decade, and motivated by a number of porous media applications, research into such coupled models has advanced modelling of processes in porous media substantially.
Building on this effort, this work presents a novel model that couples flow, heat transfer, deformation, and propagation of fractures with flow, heat transfer, and thermo-poroelasticity in the matrix. The model is based on explicit representation of fractures in the porous medium, and discretised using multi-point finite volume methods. Frictional contact and non-penetration conditions for the fractures are handled through active set methods, while a propagation criterion based on stress intensity factors governs fracture extension. Considering both forced and natural convection processes, the numerical results show the intricate nature of thermo-poromechanical fracture deformation and propagation.
△ Less
Submitted 11 December, 2020;
originally announced December 2020.
-
Hydro-mechanical simulation and analysis of induced seismicity for a hydraulic stimulation test at the Reykjanes geothermal field, Iceland
Authors:
Eirik Keilegavlen,
Laure Duboeuf,
Anna Maria Dichiarante,
Sæunn Halldórsdóttir,
Ivar Stefansson,
Marcel Naumann,
Egill Árni Guðnason,
Kristján Ágústsson,
Guðjón Helgi Eggertsson,
Volker Oye,
Inga Berre
Abstract:
The hydraulic stimulation of the well RN-34 at the Reykjanes geothermal field in Iceland caused increased seismic activity near the well. Here, we use this as a case study for investigation on how seismic analysis can be combined with physics-based simulation studies to further understand injection-induced fault reactivation. The work presents new analysis of the seismic data combined with applica…
▽ More
The hydraulic stimulation of the well RN-34 at the Reykjanes geothermal field in Iceland caused increased seismic activity near the well. Here, we use this as a case study for investigation on how seismic analysis can be combined with physics-based simulation studies to further understand injection-induced fault reactivation. The work presents new analysis of the seismic data combined with application of a recent simulation software for modeling of coupled hydromechanical processes and fault deformation caused by fluid injection. The simulation model incorporates an explicit model of the fault network based on geological characterization combined with insights from seismic analysis. The 3D faulted reservoir model is then calibrated based on injection data. Despite limited data, the work shows how seismic interpretations can be used in developing simulation models and, reciprocally, how the modeling can add to the seismic interpretations in analysis of dynamics.
△ Less
Submitted 6 November, 2020;
originally announced November 2020.
-
Fault slip in hydraulic stimulation of geothermal reservoirs: governing mechanisms and process-structure interaction
Authors:
Inga Berre,
Ivar Stefansson,
Eirik Keilegavlen
Abstract:
Hydraulic stimulation of geothermal reservoirs in low-permeability basement and crystalline igneous rock can enhance permeability by reactivation and shear dilation of existing fractures. The process is characterized by interaction between fluid flow, deformation, and the fractured structure of the formation. The flow is highly affected by the fracture network, which in turn is deformed because of…
▽ More
Hydraulic stimulation of geothermal reservoirs in low-permeability basement and crystalline igneous rock can enhance permeability by reactivation and shear dilation of existing fractures. The process is characterized by interaction between fluid flow, deformation, and the fractured structure of the formation. The flow is highly affected by the fracture network, which in turn is deformed because of hydromechanical stress changes caused by the fluid injection. This process-structure interaction is decisive for the outcome of hydraulic stimulation, and, in analysis of governing mechanisms, physics-based modeling has potential to complement field and experimental data. Here, we show how recently developed simulation technology is a valuable tool to understand governing mechanisms of hydromechanical coupled processes and the reactivation and deformation of faults. The methodology fully couples flow in faults and matrix with poroelastic matrix deformation and a contact mechanics model for the faults, including dilation because of slip. Key elements are high aspect ratios of faults and strong nonlinearities in highly coupled governing equations. Example simulations using our open-source software illustrate direct and indirect hydraulic fault reactivation and corresponding permeability enhancement. While conceptually simple, the examples illustrate the strong hydromechanical couplings and the prospects of physics-based numerical models in investigating the dynamics.
△ Less
Submitted 5 December, 2020; v1 submitted 25 August, 2020;
originally announced August 2020.
-
A fully coupled numerical model of thermo-hydro-mechanical processes and fracture contact mechanics in porous media
Authors:
Ivar Stefansson,
Inga Berre,
Eirik Keilegavlen
Abstract:
A range of phenomena in the subsurface is characterised by the interplay between coupled thermal, hydraulic and mechanical processes and deforming structures such as fractures. Modelling subsurface dynamics can provide valuable phenomenological understanding, but requires models which faithfully represent the dynamics involved; these models, therefore are themselves highly complex.
This paper pr…
▽ More
A range of phenomena in the subsurface is characterised by the interplay between coupled thermal, hydraulic and mechanical processes and deforming structures such as fractures. Modelling subsurface dynamics can provide valuable phenomenological understanding, but requires models which faithfully represent the dynamics involved; these models, therefore are themselves highly complex.
This paper presents a mixed-dimensional thermo-hydro-mechanical model designed to capture the process-structure interplay using a discrete-fracture-matrix framework. It incorporates tightly coupled thermo-hydro-mechanical processes based on laws for momentum, mass and entropy in subdomains representing the matrix and the lower-dimensional fractures and fracture intersections. The deformation of explicitly represented fractures is modelled by contact mechanics relations and a Coulomb friction law, with particular attention on coupling of fracture dilation to the governing equations in both fractures and matrix.
The model is discretised using multi-point finite volumes for the balance equations and a semismooth Newton scheme for the contact conditions and is implemented in the open source fracture simulation toolbox PorePy. Finally, simulation studies demonstrate the model's convergence, investigate process-structure coupling effects, explore different fracture dilation models and show an application of the model to a 3d geothermal pressure stimulation and long-term cooling scenario.
△ Less
Submitted 14 May, 2021; v1 submitted 14 August, 2020;
originally announced August 2020.
-
Verification benchmarks for single-phase flow in three-dimensional fractured porous media
Authors:
Inga Berre,
Wietse M. Boon,
Bernd Flemisch,
Alessio Fumagalli,
Dennis Gläser,
Eirik Keilegavlen,
Anna Scotti,
Ivar Stefansson,
Alexandru Tatomir,
Konstantin Brenner,
Samuel Burbulla,
Philippe Devloo,
Omar Duran,
Marco Favino,
Julian Hennicker,
I-Hsien Lee,
Konstantin Lipnikov,
Roland Masson,
Klaus Mosthaf,
Maria Giuseppina Chiara Nestola,
Chuen-Fa Ni,
Kirill Nikitin,
Philipp Schädle,
Daniil Svyatskiy,
Ruslan Yanbarisov
, et al. (1 additional authors not shown)
Abstract:
Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flo…
▽ More
Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flow in three-dimensional fractured porous media. The cases are specifically designed to test the methods' capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases.
△ Less
Submitted 17 February, 2020;
originally announced February 2020.
-
Finite volume discretisation of fracture deformation in thermo-poroelastic media
Authors:
Ivar Stefansson,
Inga Berre,
Eirik Keilegavlen
Abstract:
This paper presents a model where thermo-hydro-mechanical processes are coupled to a deformation model for preexisting fractures. The model is formulated within a discrete-fracture-matrix framework where the rock matrix and the fractures are considered as individual subdomains, and interaction between them takes place on the matrix-fracture interfaces. A finite volume discretisation implemented in…
▽ More
This paper presents a model where thermo-hydro-mechanical processes are coupled to a deformation model for preexisting fractures. The model is formulated within a discrete-fracture-matrix framework where the rock matrix and the fractures are considered as individual subdomains, and interaction between them takes place on the matrix-fracture interfaces. A finite volume discretisation implemented in the simulation toolbox PorePy is presented and applied in a simulation showcasing the effects of the different mechanisms on fracture deformation governed by contact mechanics, as well as their different timescales.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
PorePy: An Open-Source Software for Simulation of Multiphysics Processes in Fractured Porous Media
Authors:
Eirik Keilegavlen,
Runar Berge,
Alessio Fumagalli,
Michele Starnoni,
Ivar Stefansson,
Jhabriel Varela,
Inga Berre
Abstract:
Development of models and dedicated numerical methods for dynamics in fractured rocks is an active research field, with research moving towards increasingly advanced process couplings and complex fracture networks. The inclusion of coupled processes in simulation models is challenged by the high aspect ratio of the fractures, the complex geometry of fracture networks and the crucial impact of proc…
▽ More
Development of models and dedicated numerical methods for dynamics in fractured rocks is an active research field, with research moving towards increasingly advanced process couplings and complex fracture networks. The inclusion of coupled processes in simulation models is challenged by the high aspect ratio of the fractures, the complex geometry of fracture networks and the crucial impact of processes that completely change characteristics on the fracture-rock interface. This paper provides a general discussion of design principles for introducing fractures in simulators, and defines a framework for integrated modeling, discretization and computer implementation. The framework is implemented in the simulation software PorePy, which can serve as a flexible prototyping tool or multiphysics problems in fractured rocks. Based on a representation of the fractures and their intersections as lower-dimensional objects, we discuss data structures for mixed-dimensional meshes, formulation of multiphysics problems and discretizations that utilize existing software. We further present the implementation of these concepts in the PorePy open-source software tool, which is aimed at coupled simulation of flow and transport in three-dimensional fractured reservoirs as well as deformation of fractures and the reservoir in general. We present validation by benchmarks for flow, poroelasticity and fracture deformation in fractured porous media. The flexibility of the framework is then illustrated by simulations of fully coupled flow and transport and of injection driven deformation of fractures. All results reported herein can be reproduced by openly available simulation scripts.
△ Less
Submitted 26 August, 2019;
originally announced August 2019.
-
Call for participation: Verification benchmarks for single-phase flow in three-dimensional fractured porous media
Authors:
Inga Berre,
Wietse Boon,
Bernd Flemisch,
Alessio Fumagalli,
Dennis Gläser,
Eirik Keilegavlen,
Anna Scotti,
Ivar Stefansson,
Alexandru Tatomir
Abstract:
This call for participation proposes four benchmark tests to verify and compare numerical schemes to solve single-phase flow in fractured porous media. With this, the two-dimensional suite of benchmark tests presented by Flemisch et al. 2018 is extended to include three-dimensional problems. Moreover, transport simulations are included as a means to compare discretization methods for flow. With th…
▽ More
This call for participation proposes four benchmark tests to verify and compare numerical schemes to solve single-phase flow in fractured porous media. With this, the two-dimensional suite of benchmark tests presented by Flemisch et al. 2018 is extended to include three-dimensional problems. Moreover, transport simulations are included as a means to compare discretization methods for flow. With this publication, we invite researchers to contribute to the study by providing results to the test cases based on their applied discretization methods.
△ Less
Submitted 18 September, 2018;
originally announced September 2018.
-
Hybrid-Dimensional Finite Volume Discretizations for Fractured Porous Media
Authors:
Ivar Stefansson,
Inga Berre,
Eirik Keilegavlen
Abstract:
Over the last decade, finite volume discretizations for flow in porous media have been extended to handle situations where fractures dominate the flow. These discretizations have successfully been combined with the discrete fracture-matrix models to yield mass conservative methods capable of explicitly incorporating the impact of fractures and their geometry. When combined with a hybrid-dimensiona…
▽ More
Over the last decade, finite volume discretizations for flow in porous media have been extended to handle situations where fractures dominate the flow. These discretizations have successfully been combined with the discrete fracture-matrix models to yield mass conservative methods capable of explicitly incorporating the impact of fractures and their geometry. When combined with a hybrid-dimensional formulation, two central concerns are the restrictions arising from small cell sizes at fracture intersections and the coupling between fractures and matrix. Focusing on these aspects, we demonstrate how finite volume methods effectively can be extended to handle fractures, providing generalizations of previous work. We address the finite volume methods applying a general hierarchical formulation, facilitating implementation with extensive code reuse and providing a natural framework for coupling of different subdomains. Furthermore, we demonstrate how a Schur complement technique may be used to obtain a robust and versatile method for fracture intersection cell elimination. We investigate the accuracy of the proposed elimination method through a series of numerical simulations in 3D and 2D. The simulations, performed on fractured domains containing permeability heterogeneity and anisotropy, also demonstrate the flexibility of the hierarchical framework.
△ Less
Submitted 22 December, 2017;
originally announced December 2017.
-
Implementation of mixed-dimensional models for flow in fractured porous media
Authors:
Eirik Keilegavlen,
Alessio Fumagalli,
Runar Berge,
Ivar Stefansson
Abstract:
Models that involve coupled dynamics in a mixed-dimensional geometry are of increasing interest in several applications. Here, we describe the development of a simulation model for flow in fractured porous media, where the fractures and their intersections form a hierarchy of interacting subdomains. We discuss the implementation of a simulation framework, with an emphasis on reuse of existing disc…
▽ More
Models that involve coupled dynamics in a mixed-dimensional geometry are of increasing interest in several applications. Here, we describe the development of a simulation model for flow in fractured porous media, where the fractures and their intersections form a hierarchy of interacting subdomains. We discuss the implementation of a simulation framework, with an emphasis on reuse of existing discretization tools for mono-dimensional problems. The key ingredients are the representation of the mixed-dimensional geometry as a graph, which allows for convenient discretization and data storage, and a non-intrusive coupling of dimensions via boundary conditions and source terms. This approach is applicable for a wide class of mixed-dimensional problems. We show simulation results for a flow problem in a three-dimensional fracture geometry, applying both finite volume and virtual finite element discretizations.
△ Less
Submitted 20 December, 2017;
originally announced December 2017.
-
PorePy: An Open-Source Simulation Tool for Flow and Transport in Deformable Fractured Rocks
Authors:
Eirik Keilegavlen,
Alessio Fumagalli,
Runar Berge,
Ivar Stefansson,
Inga Berre
Abstract:
Fractures are ubiquitous in the subsurface and strongly affect flow and deformation. The physical shape of the fractures, they are long and thin objects, puts strong limitations on how the effect of this dynamics can be incorporated into standard reservoir simulation tools. This paper reports the development of an open-source software framework, termed PorePy, which is aimed at simulation of flow…
▽ More
Fractures are ubiquitous in the subsurface and strongly affect flow and deformation. The physical shape of the fractures, they are long and thin objects, puts strong limitations on how the effect of this dynamics can be incorporated into standard reservoir simulation tools. This paper reports the development of an open-source software framework, termed PorePy, which is aimed at simulation of flow and transport in three-dimensional fractured reservoirs, as well as deformation of the reservoir due to shearing along fracture and fault planes. Starting from a description of fractures as polygons embedded in a 3D domain, PorePy provides semi-automatic gridding to construct a discrete-fracture-matrix model, which forms the basis for subsequent simulations. PorePy allows for flow and transport in all lower-dimensional objects, including planes (2D) representing fractures, and lines (1D) and points (0D), representing fracture intersections. Interaction between processes in neighboring domains of different dimension is implemented as a sequence of couplings of objects one dimension apart. This readily allows for handling of complex fracture geometries compared to capabilities of existing software. In addition to flow and transport, PorePy provides models for rock mechanics, poro-elasticity and coupling with fracture deformation models. The software is fully open, and can serve as a framework for transparency and reproducibility of simulations. We describe the design principles of PorePy from a user perspective, with focus on possibilities within gridding, covered physical processes and available discretizations. The power of the framework is illustrated with two sets of simulations; involving respectively coupled flow and transport in a fractured porous medium, and low-pressure stimulation of a geothermal reservoir.
△ Less
Submitted 21 December, 2017; v1 submitted 1 December, 2017;
originally announced December 2017.
-
Benchmarks for single-phase flow in fractured porous media
Authors:
Bernd Flemisch,
Inga Berre,
Wietse Boon,
Alessio Fumagalli,
Nicolas Schwenck,
Anna Scotti,
Ivar Stefansson,
Alexandru Tatomir
Abstract:
This paper presents several test cases intended to be benchmarks for numerical schemes for single-phase fluid flow in fractured porous media. A number of solution strategies are compared, including a vertex and a cell-centered finite volume method, a non-conforming embedded discrete fracture model, a primal and a dual extended finite element formulation, and a mortar discrete fracture model. The p…
▽ More
This paper presents several test cases intended to be benchmarks for numerical schemes for single-phase fluid flow in fractured porous media. A number of solution strategies are compared, including a vertex and a cell-centered finite volume method, a non-conforming embedded discrete fracture model, a primal and a dual extended finite element formulation, and a mortar discrete fracture model. The proposed benchmarks test the schemes by increasing the difficulties in terms of network geometry, e.g. intersecting fractures, and physical parameters, e.g. low and high fracture-matrix permeability ratio as well as heterogeneous fracture permeabilities. For each problem, the results presented by the participants are the number of unknowns, the approximation errors in the porous matrix and in the fractures with respect to a reference solution, and the sparsity and condition number of the discretized linear system. All data and meshes used in this study are publicly available for further comparisons.
△ Less
Submitted 5 January, 2017;
originally announced January 2017.