Researcher profile

Eirik Keilegavlen

Eirik Keilegavlen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2023arXiv

PoroTwin: A digital twin for a FluidFlower rig

We present a framework for integrated experiments and simulations of tracer transport in heterogeneous porous media using digital twin technology. The physical asset in our setup is a meter-scale FluidFlower rig. The digital twin consists of a traditional physics-based forward simulation tool and a correction technique which compensates for mismatches between simulation results and observations. The latter augments the range of the physics-based simulation and allows us to bridge the gap between simulation and experiments in a quantitative sense. We describe the setup of the physical and digital twin, including data transfer protocols using cloud technology. The accuracy of the digital twin is demonstrated on a case with artificially high diffusion that must be compensated by the correction approach, as well as by simulations in geologically complex media. The digital twin is then applied to control tracer transport by manipulating fluid injection and production in the experimental rig, thereby enabling two-way coupling between the physical and digital twins.

preprint2022arXiv

A posteriori error estimates for hierarchical mixed-dimensional elliptic equations

Mixed-dimensional elliptic equations exhibiting a hierarchical structure are commonly used to model problems with high aspect ratio inclusions, such as flow in fractured porous media. We derive general abstract estimates based on the theory of functional a posteriori error estimates, for which guaranteed upper bounds for the primal and dual variables and two-sided bounds for the primal-dual pair are obtained. We improve on the abstract results obtained with the functional approach by proposing four different ways of estimating the residual errors based on the extent the approximate solution has conservation properties, i.e.: (1) no conservation, (2) subdomain conservation, (3) grid-level conservation, and (4) exact conservation. This treatment results in sharper and fully computable estimates when mass is conserved either at the grid level or exactly, with a comparable structure to those obtained from grid-based a posteriori techniques. We demonstrate the practical effectiveness of our theoretical results through numerical experiments using four different discretization methods for synthetic problems and applications based on benchmarks of flow in fractured porous media.

preprint2022arXiv

Impact of deformation bands on fault-related fluid flow in field-scale simulations

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.

preprint2021arXiv

Modelling and discretization of flow in porous media with thin, full-tensor permeability inclusions

When modelling fluid flow in fractured reservoirs, it is common to represent the fracturesas lower-dimensional inclusions embedded in the host medium. Existing discretizationsof flow in porous media with thin inclusions assume that the principal directions of theinclusion permeability tensor are aligned with the inclusion orientation. While this mod-elling assumption works well with tensile fractures, it may fail in the context of faults,where the damage zone surrounding the main slip surface may introduce anisotropy thatis not aligned with the main fault orientation. In this paper, we introduce a generalizeddimensional reduced model which preserves full-tensor permeability effects also in theout-of-plane direction of the inclusion. The governing equations of flow for the lower-dimensional objects are obtained through vertical averaging. We present a framework fordiscretization of the resulting mixed-dimensional problem, aimed at easy adaptation ofexisting simulation tools. We give numerical examples that show the failure of existingformulations when applied to anisotropic faulted porous media, and go on to show theconvergence of our method in both 2D and 3D

preprint2020arXiv

An Introduction to Multi-Point Flux (MPFA) and Stress (MPSA) Finite Volume Methods for Thermo-Poroelasticity

We give a unified introduction to the MPFA- and MPSA-type finite volume methods for Darcy flow and poro-elasticity, applicable to general polyhedral grids. This leads to a more systematic perspective of these methods than has been exposed in previous texts, and we therefore refer to this discretization family as the MPxA methods. We apply this MPxA framework to also define a consistent finite-volume discretization of thermo-poro-elasticity. In order to make the exposition accessible to a wide audience, we avoid much of the technical notation which is used in the research literature, and compensate for this by an expanded summary and literature review of the main properties of the MPxA methods. We close the chapter by a section containing applications to problems with complex geometries and non-linear physics.

preprint2020arXiv

Finite volume discretisation of fracture deformation in thermo-poroelastic media

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.

preprint2020arXiv

Numerical modeling of wing crack propagation accounting for fracture contact mechanics

As a consequence of shearing, wing cracks can emerge from pre-existing fractures. The process involves the interaction of sliding of the existing fracture surfaces and the tensile material failure that creates wing cracks. This work devises a numerical model to investigate how wing cracks emerge, propagate and connect pre-existing fractures under shear processes. A mathematical and numerical model for wing crack propagation based on linear elastic fracture mechanics that also accounts for fracture contact mechanics is presented. Computational efficiency is ensured by an adaptive remeshing technique. The numerical model is verified and validated through a comparison of the analytical and experimental results. Additional numerical examples illustrate the performance of the method for complex test cases where wing-cracks develop for multiple pre-existing and interacting fractures.