Numerical Modeling of the Effect of Carbon Dioxide Sequestration on the Rate of Pressure Solution Creep in Limestone: Preliminary Results

Abstract

When carbon dioxide (CO2 ) is injected into an aquifer or a depleted geological reservoir, its dissolution into solution results in acidification of the pore waters. As a consequence, the pore waters become more reactive, which leads to enhanced dissolution-precipitation processes and a modification of the mechanical and hydrological properties of the rock. This effect is especially important for limestones given that the solubility and reactivity of carbonates is strongly dependent on pH and the partial pressure of CO2 . The main mechanism that couples dissolution, precipitation and rock matrix deformation is commonly referred to as intergranular pressure solution creep (IPS) or pervasive pressure solution creep (PSC). This process involves dissolution at intergranular grain contacts subject to elevated stress, diffusion of dissolved material in an intergranular fluid, and precipitation in pore spaces subject to lower stress. This leads to an overall and pervasive reduction in porosity due to both grain indentation and precipitation in pore spaces. The percolation of CO2 -rich fluids may influence on-going compaction due to pressure solution and can therefore potentially affect the reservoir and its long-term CO2 storage capacity. We aim at quantifying this effect by using a 2D numerical model to study the coupling between dissolution-precipitation processes, local mass transfer, and deformation of the rock over long time scales. We show that high partial pressures of dissolved CO2 (up to 30 MPa) significantly increase the rates of compaction by a factor of ~50 to ~75, and also result in a concomitant decrease in the viscosity of the rock matrix.