TY - JOUR

T1 - Computational topology-based characterization of pore space changes due to chemical dissolution of rocks

AU - Lisitsa, Vadim

AU - Bazaikin, Yaroslav

AU - Khachkova, Tatyana

N1 - Publisher Copyright: © 2020 Elsevier Inc. Copyright: Copyright 2020 Elsevier B.V., All rights reserved.

PY - 2020/12/1

Y1 - 2020/12/1

N2 - In this paper, we present an algorithm for the numerical simulation of reactive transport at the pore scale to facilitate observation of pore space and rock matrix evolution. Moreover, simulation at the pore scale opens up the possibility of estimating changes in the transport properties of rocks, such as permeability and tortuosity. To quantitatively analyze pore space evolution, we developed a numerical algorithm that can be used to construct persistence diagrams of the connectivity components for pore space and the rock matrix, which characterize the topology evolution during rock matrix dissolution. Introducing the “bottle-neck” metric in the space of the persistence diagrams, we cluster the numerical experiments in terms of similarities in topology evolution. We demonstrate that the application of this metric to the persistence diagrams allowed us to distinguish topologically different dissolution scenarios, for instance, the formation of a dissolution front near the inlet, homogeneous dissolution of the matrix inside the core sample, and formation of wormholes. We illustrate that the differences in topology evolution lead to cross-correlations among the transport properties of rocks (porosity-permeability-tortuosity).

AB - In this paper, we present an algorithm for the numerical simulation of reactive transport at the pore scale to facilitate observation of pore space and rock matrix evolution. Moreover, simulation at the pore scale opens up the possibility of estimating changes in the transport properties of rocks, such as permeability and tortuosity. To quantitatively analyze pore space evolution, we developed a numerical algorithm that can be used to construct persistence diagrams of the connectivity components for pore space and the rock matrix, which characterize the topology evolution during rock matrix dissolution. Introducing the “bottle-neck” metric in the space of the persistence diagrams, we cluster the numerical experiments in terms of similarities in topology evolution. We demonstrate that the application of this metric to the persistence diagrams allowed us to distinguish topologically different dissolution scenarios, for instance, the formation of a dissolution front near the inlet, homogeneous dissolution of the matrix inside the core sample, and formation of wormholes. We illustrate that the differences in topology evolution lead to cross-correlations among the transport properties of rocks (porosity-permeability-tortuosity).

KW - Betti numbers

KW - Computational topology

KW - Finite-differences

KW - Level-set

KW - Persistence diagram,

KW - Rock dissolution

KW - Persistence diagram

KW - POROUS-MEDIA

KW - MECHANISMS

KW - SCALE

KW - LEVEL-SET METHODS

KW - FLOW

UR - http://www.scopus.com/inward/record.url?scp=85087591561&partnerID=8YFLogxK

U2 - 10.1016/j.apm.2020.06.037

DO - 10.1016/j.apm.2020.06.037

M3 - Article

AN - SCOPUS:85087591561

VL - 88

SP - 21

EP - 37

JO - Applied Mathematical Modelling

JF - Applied Mathematical Modelling

SN - 0307-904X

ER -