Carbon dioxide (CO2) injection into deep saline aquifers is governed by a number of physico-chemical processes including mineral dissolution and precipitation, multiphase fluid flow, and capillary trapping. These processes can be coupled, however, the impact of fluid-rock reaction on the multiphase flow properties is difficult to study and is not simply correlated to variation in rock porosity. We observed the impact of rock mineral dissolution on multiphase flow properties in two carbonate rocks with distinct pore structures. The Ketton carbonate was an ooidal limestone with a distinct bimodal pore structure whereas the Estaillades limestone was a bioclastic limestone with a wide range of pore sizes. Observations of steady state N2-water relative permeability and residual trapping were obtained at 100 bars fluid pressure and 22°C, with X-ray tomography used to estimate fluid saturation. These tests alternated with steps in which mineral was uniformly dissolved into solution from the rock cores using an aqueous solution with a temperature controlled acid. Eight alternating sequences of dissolution and flow measurement were performed, with on average 0.5% of the mass of the rocks dissolved at each stage. A sequence of mercury injection capillary pressure measurements were conducted on a parallel set of samples undergoing the same treatment to characterize the evolving pore size distribution and corresponding capillary pressure characteristics. Variations in the multiphase flow properties were observed to correspond to the changes in the underlying pore structure. In the Ketton carbonate, dissolution resulted in an increase of the fraction of pore volume made up by the smallest pores and a corresponding increase in the fraction made up by the largest pores. This resulted in a systematic increase in the relative permeability to the nonwetting phase and decrease in relative permeability of the wetting phase. There was also a modest but systematic decrease in residual trapping. In the Estaillades carbonate, dissolution resulted in an increase in the fraction of pore volume made up by pores in the central range of the initial pore size distribution, and a corresponding decrease in the fraction made up by both the smallest and largest pores. This resulted in a decrease in the relative permeability to both the wetting and nonwetting fluid phases and no discernible impact on the residual trapping. In summary, the impact of rock matrix dissolution will be strongly dependent on the impact of that dissolution on the underlying pore structure of the rock. However, if the variation in pore structure can be observed or estimated with modelling, then it should be possible to estimate the impacts on multiphase flow properties.
