Microtomographic characterization of dissolution-induced local porosity changes including fines migration in carbonate rock

The microstructure of carbonate rocks experiences substantial changes under reactive processes, in particular chemical dissolution and deposition, including dissolution-released-fines migration occurring during acidizing. A better understanding of such changes at the pore scale and their influences on rock properties is of great value for the effective design and implementation of reactive processes in carbonate reservoirs. In this work, we demonstrate the use of X-ray micro-computed tomography (micro-CT) to quantitatively investigate the local porosity changes in a meso-/microporous carbonate core sample during chemical dissolution. A reactive flooding experiment in a core sample by a nonacidic solution is designed such that changes in pore space from before to after the reactant injection could be imaged in exactly the same locations with micro-CT at a resolution of less than 5 μm. A methodology with three-phase segmentation and 2D histograms of image intensity is used to quantify distributions of the evolution of each image voxel. This technique allows the incorporation of microporosity into the calculation of the evolution regions, including the migration of fines, to accurately quantify the evolution scenarios. The micro-CT images reveal a quasiuniform dissolution pattern and allow characterizing the accompanying migration of fines within the core sample. The 3D pore networks are derived from the image data, which quantify changes in network structure and the pore geometry. The 2D histograms of image intensity derived from the pre- and post-dissolution images show quantitatively how macro- and micropores are enlarged by dissolution close to the inlet, whereas the deposition of fines mainly occurs in pores far from the inlet boundary. These results can explain why permeability of the sample initially decreases and then increases when injection time increases. Pore-surface area between each region is computed on the basis of the spatially resolved voxel evolution scenarios. This allows calculation of local distribution of reactive surface area, which, in turn, will assist in the prediction of local reaction rates in reactive flow simulators.