Explore the vast range of titles available.

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 open-source simulation software PorePy, which can serve as a flexible prototyping tool for 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 grids, formulation of multiphysics problems, and discretizations that utilize existing software. We further present a Python 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 porous media. The flexibility of the framework is then illustrated by simulations of non-linearly coupled flow and transport and of injection-driven deformation of fractures. All results can be reproduced by openly available simulation scripts.

The effect of heterogeneity induced by highly permeable fracture networks on viscous miscible fingering in porous media is examined using high-resolution numerical simulations. We consider the planar injection of a less viscous fluid into a two-dimensional fractured porous medium that is saturated with a more viscous fluid. This problem contains two sets of fundamentally different preferential flow regimes; the first is caused by the viscous fingering, and the second is due to the permeability contrasts between the fractures and the rock matrix. We study the transition from the regime where the flow is dominated by the viscous instabilities, to the regime where the heterogeneity induced by the fractures define the flow paths. Our findings reveal that even minor permeability differences between the rock matrix and fractures significantly influence the behavior of viscous fingering. The interplay between the viscosity contrast and permeability contrast leads to the preferential channeling of the less viscous fluid through the fractures. Consequently, this channeling process stabilizes the displacement front within the rock matrix, ultimately suppressing the occurrence of viscous fingering, particularly for higher permeability contrasts. We explore three fracture geometries: two structured and one random configuration and identify a complex interaction between these geometries and the development of unstable flow. While we find that the most important factor determining the effect of the fracture network is the ratio of fluid volume flowing through the fractures and the rock matrix, the exact point for the cross-over regime is dependent on the geometry of the fracture network.

We present a novel mixed-dimensional method for generating unstructured polyhedral grids that conform to prescribed geometric objects in arbitrary dimensions. Two types of conformity are introduced: (i) control-point alignment of cell centroids to accurately represent horizontal and multilateral wells or create volumetric representations of fracture networks, and ii) boundary alignment of cell faces to accurately preserve lower dimensional geological objects such as layers, fractures, faults, and/or pinchouts. The prescribed objects are in this case assumed to be lower dimensional, and we create a grid hierarchy in which each lower dimensional object is associated with a lower dimensional grid. Further, the intersection of two objects is associated with a grid one dimension lower than the objects. Each grid is generated as a clipped Voronoi diagram, also called a perpendicular bisector (PEBI) grid, for a carefully chosen set of generating points. Moreover, each grid is generated in such a way that the cell faces of a higher dimensional grid conform to the cells of all lower dimensional grids. We also introduce a sufficient and necessary condition which makes it easy to check if the sites for a given perpendicular bisector grid will conform to the set of prescribed geometric objects.

Subsurface storage of CO2 is predicted to rise exponentially in response to the increasing levels of CO2 in the atmosphere. Large-scale CO2 injections into the subsurface require understanding of the potential for fluid flow through faults to mitigate risk of leakage. Here, we study how to obtain effective permeability of deformation bands in the damage zone of faults. Deformation bands are relatively small, low permeability features that can have a significant effect on flow dynamics, however, the discrepancy of scales is a challenge for field-scale simulation. A new analytical upscaling model is proposed in order to overcome some of the shortcomings of conventional upscaling approaches for heterogeneous porous media. The new model captures the fine-scale impact of deformation bands on fluid flow in the near-fault region, and can be derived from knowledge of large-scale fault properties. To test the accuracy of the model it is compared to fine-scale numerical simulations that explicitly include individual deformation bands. For a wide range of different stochastically generated deformation bands networks, the upscaling model shows improved estimate of effective permeability compared to conventional upscaling approaches. By applying the upscaling model to a full-field simulation of the Smeaheia storage site in the North Sea, we show that deformation bands with a permeability contrast higher than three orders of magnitude may act as an extra layer of protection from fluid flow through faults.

Fluid-induced slip of fractures is characterized by strong multiphysics couplings. Three physical processes are considered: Flow, rock deformation and fracture deformation. The fractures are represented as lower-dimensional objects embedded in a three-dimensional domain. Fluid is modeled as slightly compressible, and flow in both fractures and matrix is accounted for. The deformation of rock is inherently different from the deformation of fractures; thus, two different models are needed to describe the mechanical deformation of the rock. The medium surrounding the fractures is modeled as a linear elastic material, while the slip of fractures is modeled as a contact problem, governed by a static-dynamic friction model. We present an iterative scheme for solving the non-linear set of equations that arise from the models, and suggest how the step parameter in this scheme should depend on the shear modulus and mesh size.

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.

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.

The effect of heterogeneities induced by highly permeable fracture networks on viscous miscible fingering in porous media is examined using high-resolution numerical simulations. We consider the planar injection of a less viscous fluid into a two-dimensional fractured porous medium which is saturated with a more viscous fluid. This problem contains two sets of fundamentally different preferential flow regimes; the first is caused by the viscous fingering and the second is due to the permeability contrasts between the fractures and the rock matrix. We study the transition from the regime where the flow is dominated by the viscous instabilities, to the regime where the heterogeneities induced by the fractures define the flow paths. Our findings reveal that even minor permeability differences between the rock matrix and fractures significantly influence the behavior of viscous fingering. The interplay between the viscosity contrast and permeability contrast leads to the preferential channeling of the less viscous fluid through the fractures. Consequently, this channeling process stabilizes the displacement front within the rock matrix, ultimately suppressing the occurrence of viscous fingering, particularly for higher permeability contrasts. We explore three fracture geometries; two structured and one random configuration, and identify a complex interaction between these geometries and the development of unstable flow. While we find that the most important factor determining the effect of the fracture network is the ratio of fluid volume flowing through the fractures and the rock matrix, the exact point for the cross-over regime is dependent on the geometry of the fracture network.

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.

A fractured poroelastic body is considered where the opening of the fractures is governed by a nonpenetration law while slip is described by a Coulomb-type friction law. This physical model results in a nonlinear variational inequality problem. The variational inequality is rewritten as a complimentary function, and a semismooth Newton method is used to solve the system of equations. For the discretization, we use a hybrid scheme where the displacements are given in terms of degrees of freedom per element, and an additional Lagrange multiplier representing the traction is added on the fracture faces. The novelty of our method comes from combining the Lagrange multiplier from the hybrid scheme with a finite volume discretization of the poroelastic Biot equation, which allows us to directly impose the inequality constraints on each subface. The convergence of the method is studied for several challenging geometries in 2d and 3d, showing that the convergence rates of the finite volume scheme do not deteriorate when it is coupled to the Lagrange multipliers. Our method is especially attractive for the poroelastic problem because it allows for a straightforward coupling between the matrix deformation, contact conditions, and fluid pressure.