Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis- part 1: Assessment of pore-scale mineral dissolution and deposition.

The application of X-ray micro-computed tomography (μ-CT) for quantitatively characterizing reactive-flow induced pore structure evolution including local particle detachment, displacement and deposition in carbonate rocks is investigated. In the studies conducted in this field of research, the experimental procedure has involved alternating steps of imaging and ex-situ core sample alteration. Practically, it is impossible to return the sample, with micron precision, to the same position and orientation. Furthermore, successive images of a sample in pre- and post-alteration states are usually taken at different conditions such as different scales, resolutions and signal-to-noise ratios. These conditions accompanying with subresolution features in the images make voxel-by-voxel comparisons of successive images problematic. In this paper, we first address the respective challenges in voxel-wise interpretation of successive images of carbonate rocks subject to reactive flow. Reactive coreflood in two carbonate cores with different rock types are considered. For the first rock, we used the experimental and imaging results published by Qajar et al. (2013) which showed a quasi-uniform dissolution regime. A similar reactive core flood was conducted in the second rock which resulted in wormhole-like dissolution regime. We particularly examine the major image processing operations such as transformation of images to the same grey-scale, noise filtering and segmentation thresholding and propose quantitative methods to evaluate the effectiveness of these operations in voxel-wise analysis of successive images of a sample. In the second part, we generalize the methodology based on the three-phase segmentation of normalized images, microporosity assignment and 2D histogram of image intensities to estimate grey-scale changes of individual image voxels for a general case where the greyscale images are segmented into arbitrary number of phases. The results show that local (voxel-based) porosity changes can be decomposed into local mineral dissolution and deposition. Moreover, it is found that the microporosity evolutions are differently distributed in the samples after the reactive coreflood experiments. In the last part of the paper, for the case of quasi-uniform dissolution, we combine the tomographic images with numerical calculations of permeability along the core to characterize the relationship between changes in permeability and the fractions of the mineral dissolved and deposited. A consistency is found between the calculated longitudinal permeability changes and the quantified distribution of mineral dissolved and deposited along the sample.