Abstract 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).

中文翻译：

---

基于计算拓扑的岩石化学溶解引起的孔隙空间变化表征

摘要 在本文中，我们提出了一种孔隙尺度反应输运数值模拟算法，以方便观察孔隙空间和岩石基质演化。此外，孔隙尺度的模拟开辟了估计岩石传输特性变化的可能性，例如渗透率和曲折度。为了定量分析孔隙空间演化，我们开发了一种数值算法，该算法可用于构建孔隙空间和岩石基质连通性成分的持续图，表征岩石基质溶解过程中的拓扑演化。在持久图的空间中引入“瓶颈”度量，我们根据拓扑演化的相似性对数值实验进行聚类。我们证明，将该度量应用于持久性图使我们能够区分拓扑不同的溶解情况，例如，入口附近溶解前沿的形成、岩心样品内基质的均匀溶解以及虫洞的形成。我们说明拓扑演化的差异导致岩石输运性质（孔隙度-渗透率-曲折度）之间的互相关。