Dynamic modelling of overprinted low-permeability fault cores and surrounding damage zones as lower dimensional interfaces for multiphysics simulations

• Published in: Computers & geosciences Vol. 150; p. 104719
• Main Authors: Poulet, Thomas, Lesueur, Martin, Kelka, Ulrich
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.05.2021
• Subjects: Dynamic mixed permeability; Fault overprinting; Flow in fractured porous media; Lower dimensional interface; Flow in fractured porous media; Fault overprinting; Lower dimensional interface; Dynamic mixed permeability
• ISSN: 0098-3004; 1873-7803
• DOI: 10.1016/j.cageo.2021.104719

In the modelling of subsurface fluid flow, faults are dominant features since they can act as fluid pathways, baffles or barriers. Special emphasis is therefore placed in representing them in a numerically efficient manner and the use of lower dimensional domains has become prevalent to simulate higher permeability features like fractures. Such features, however, only represent some of the components of natural fault networks, which can also include rather impermeable fault gouges surrounded by higher permeability damage zones for instance. Here we present a numerical approach to simulate such systems at a large scale, where the thickness of those features makes it advantageous to represent them as discrete rather than continuous domains, using lower dimensional interfaces in a conforming mesh. Benchmarks show excellent agreement with equivalent continuous simulations, regardless of the fault thickness or permeability, both for flow conduits and baffles. This approach can also account for the overprinting of faults with different permeabilities, as well as their dynamic evolution, which we illustrate with an example of trap charging. This work demonstrates the applicability of the approach to simulate fluid flow in faulted environments of various permeabilities and we discuss how those results can easily be extended to account for multi-physical processes. •Modelling damage zones around low-permeability fault core as zero thickness elements.•Simulating fluid flow in fault zones transitioning between baffle and conduit.•Accounting for overprinting of faults with different permeabilities.